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1.0 


INFRODUCTION AND SUMMARY 


This report documents the analyses and investigations performed 
by Axiomatix under Contract NAS 9-15387C for NASA Johnson Space Center. 
This report and its companion report performed under Contract NAS 9-15387B 
present the results of work performed to analyze areas of a GPS navigation 
system to be implemented on the Shuttle that would have the greatest 
impact on performance. Also, this report presents results of the inves- 
tigation of the integration problems. The various aspects of GPS receiver 
design trade-offs that were evaluated as part of the analyses and inves- 
tigations are included in this report. 

The analysis of the Shuttle/GPS link is of foremost importance. 

It is necessary to demonstrate a sufficiently high signal-to-noise ratio 
for this link at this pont in time so that a high degree of confidence 
can be established for meeting performance requirements when the flight 
system is implemented. The approximately 3 dB of minimum margin found 
in the link analysis satisfies this goal. 

Since the Shuttle GPS antennas must be located remotely from the 
receiver, a preamplifier is necessary. This preamplifier determines the 
system noise temperature and thus strongly influences the link margin. 
Consequently, trade-offs and performance for Shuttle GPS application 
has bren investigated. Furthermore, performance of the system as influ- 
enced by redundancy switching considerations has been analyzed. 

A discussion of several GPS receiver architecture trade-offs is 
presented. The discussion is centered on the sequential receiver, which 
is recommended as the baseline receiver architecture. Three receiver 
frequency conversion plans are discussed and double conversion is high- 
lighted as an attractive choice. 

The Shuttle GPS system must operate in a complex RF environment, 
in which several on-board RF emitters may operate simultaneously. Con- 
sequently, and RFI/EMI analysis has been performed which evaluates the 
Shuttle RF harmonics and intermods that fall within the GPS receiver 
bandwidth. The analysis shows that no serious RFI/EMI problems exist. 

An analysis of TACAN/GPS EMC has shown that 28 dB of isolation between 
the TACAN antenna and the GPS antenna is desirable. It is pointed out 
that if the requirement for ionospheric delay correction is relaxed, the 



isolation can be smaller since only the L2 frequency is affected. Opera- 
t1o«'a1 procedures can also relax the isolation requirement. 

The GPS PN code acquisition has been analyzed extensively. This 
analysis has used the results of the link budget analysis to optimize 
the type of acquisition technique (sequential} and to optimize the receixe. 
parameters such as bandwidth and threshold that affect the acquisition 
time. 

The GPS signal carrier must be acquired for data detection and 
pseudo range rate measurements. In the low GPS signal level environment 
of the Shuttle, it is important that the receiver processing be optimized 
for fast carrier acquisition. A technique for potentially improving the 
carrier acquisition performance for the Shuttle GPS is presented and 
analyzed in detail. An improvement over conventional coherent receiver 
techniques is shown. 

Since the receiver clock strongly affects both GPS carrier and 
code acquisition performance, a clock model has been developed and analyzed. 
This model is useful in specifying the performance for the Shuttle GPS 
receiver clock. 

The Shuttle GPS system performance will, to a large extent, 
depend on elements of the overall GPS system beyond the control of NASA. 

For this reason, an interface control document (ICO) between NASA and 
the GPS Joint Program Office is necessary. A preliminary outline of this 
ICD is presented. 
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2.0 GPS/SHUTTLE LINK BUDGETS 

Summary 

The link budgets for the GPS satellite to the Space Shuttle 
(nominal orbit) have been calculated. Tables 1 through 12 are the 
detailed budgets for ranging, carrier tracking and data detection for 
LI and L2 frequencies and P and C/A codes. In general, the parameter val- 
ues are either worst-case specified values, best engineering predictions 
based on breadboard measurements, or survey of available equipment or com- 
ponents. The system effective noise temperature is dircussed in detail 
below but, since it is so strongly influenced by preamp performance. 
Section 3.0 has been devoted to preamplifier design and analysis. 

The link margins are summarized in Table 13. In general, the 
in-atmosphere performance should be essentially the same since the 
SS-6PS-2000 System Specification for the NAVSTAR Global Positioning Sys- 
tem specifies LI and L2 signal levels at the earth's surface of -163 and 
-166 dBW for a 0 dB circularly polarized antenna. These signal levels 
are within 0.4 dB of the GPS/Shuttle link for a -4 dB antenna gain. 

An explanation of some of the budget parameter values follows below. 

Satellite EIRP 

These values are for the edge of earth coverage as specified in 
the GPS Space Segment Specification CID-SV-IOH. The actual performance 
of the two on-orbit spacecraft has been reported to be better than the 
specification values. 

Space Loss 

Two space loss numbers of -184.6 and -182,4 dB for Ll and L2, 
on the orbital geometry shown in Figure 1. 
of sight (LOS) is given by; 

• 26.659 cos |s1n-' 

= 25,795 km 

by L = 92.45 + 20 log f (GHz) + log R (km) 

= 92.45 + 20 log 1.575 + 20 log 25,795 
= 184.6 dB. 


respectively, are based 
From Figure 1, the line 

R = 26,650 cos o 
The space loss is given 
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Table 1. Baseline Link Budget for Range Measurement, Ll-P 


Link Ll-P, Range 
Upper Antenna 

Orbital Geometry: Path Tangential to Orbit 


Parameter 

VaJ 

lue 


Explanation 

Satellite EIRP 

23.8 

dr- 

EOE, 

Spec. CID-SV-lOH 

Space Loss 

-184.6 

dB 

Path 

Tangent to Orbit 

Pointin: Loss 

- 0.4 

dB 



Polarization Loss 

- 0.4 

dB 



Atmospheric Loss 

0 




Shuttle Antenna Gain 

- 4.0 

dB 



Received Power 

-165.6 

dBW 

At Antenna Terminals 

System Noise Temperature 

28.5 

dBK 

^ C c 

* 705®K (at antenna 



sys 

terminals) 

Boltzmann's Constant 

-228.6 

dB-W/K/Hz 



Noise Spectral Density 

-200.1 

dB-W/Hz 



C/Nq 

34.5 

dB-Hz 



Required C/N^ (Theoretical) 

29.6 

dB-Hz 

®Range ' ^ 

Implementation Loss 

2.0 

dB 



Required C/Nq 

31.6 

dB-Hz 



Link Margin 

2.9 

dB 
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Table 2. Baseline Link Budget for Range Measurement, LI -C/A 


Link LI -C/A, Range 
Upper Antenna 

Orbital Geometry: Path Tangential to Orbit 


Parameter 

Value 

Explanation 

Satellite EIRP 

26.8 dBW 

EOE, Spec. CIO-SV-lOH 

Space Loss 

-184.6 dB 

Path Tangent to Orbit 

Pointing Loss 

- 0.4 dB 


Polarization Loss 

- 0.4 dB 


Atmospheric Loss 

0 


Shuttle Ai.tenna Gain 

- 4.0 dB 


Received Power 

-162.6 dBW 

At Antenna Terminals 

System Noise Temperature 
Boltzmann's Constant 

28.5 dBK 
-228.6 dB-W/K/Hz 

T = 705°K (at antenna 

^ terminals) 

Noise Spectral Density 

-200.1 dB-W/Hz 


C/Nq 

37.5 dB-Hz 


Required C/Nq (Theoretical) 

29.6 dB-Hz 

°Range ~ 

Implementation Loss 

2.0 dB 


Required C/Nq 

31.6 dB-Hz 



Link Margin 


5.9 dB 



i 


Table 3. Baseline Link Budget for Range Measurement, L2-P 


Link L2-P, Range 
Upper Antenna 

Orbital Geometry: Path Tangential to Orbit 


Parameter 

Value 

Explanation 

Satellite EIRP 

19.1 dBW 

EOE, Spec. CID-SV-lOH 

Space Loss 

-182.4 dB 

Path Tangent to Orbit 

Pointing Loss 

- 0.4 dB 


Polarization Loss 

- 0.4 dB 


Atmospheric Loss 

0 


Shuttle Antenna Gain 

- 4.0 dB 


Received Power 

-168.1 dBW 

At Antenna Terminals 

System Noise Temperature 
Boltzmann's Constant 

28.5 dBK 

-228.6 dB-W/K/Hz 

T = 705“K (at antenna 

^ terminals) 

Noise Spectral Density 

-200.1 dB/W/Hz 


C/Nq 

32.0 dB-Hz 


Required C/Nq (Theoretical) 

27.5 dB-Hz 

'^Range * ^ 

Implementation Loss 

2.0 dB 


Required C/Nq 

29.5 dB 


Link Margin 

2.5 dB 
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Table 4. Baseline Link Budget for Range Measurement, L2-C/A 


Link L2-C/A, Range 
Upper Antenna 

Orbital Geometry: Path Tangential to Orbit 


Parameter 

Value 

Explanation 

Satellite EIRP 

19.1 dBW 

EOE, Spec. CID-SV-lOH 

Space Loss 

-182.4 dB 

Path Tangent to Orbit 

Pointing Loss 

- 0.4 dB 


Polarization Loss 

- 0.4 dB 


Atmospheric Loss 

0 


Shuttle Antenna Gain 

- 4.0 dB 


Received Power 

-168.1 dBW 

At Antenna Terminals 

System Noise Temperature 

28.5 dBK 

T = TOB^K (at antenna 

terminals) 

Boltzmann's Constant 

-228.6 dB-W/K/Hz 


Noise Spectral Density 

-200.1 dB-W/Hz 


C/Nq 

32.0 dB-Hz 


Required C/Nq (Theoretical) 

27.5 dB-Hz 

'■’Range ~ 

Implementation Loss 

2.0 dB 


Required C/Nq 

29.5 dB-Hz 


Link Margin 

2.5 dB 




1 ■ 


8 


Table 5. Baseline Link Budget for Ca>^rier Tracking, Ll-P 


Link Ll-P, Carrier 
Upper Antenna 

Orbital Geometry: Path Tangential to Orbit 


Parameter 

Value 

Explanation 

Satellite EIRP 

23.8 dBW 

EOE, Spec. CID-SV-lOH 

Space Loss 

-184.6 dB 

Path Tangent to Orbit 

Pointing Loss 

- 0.4 dB 


Polarization Loss 

- 0.4 dB 


Atmospheric Loss 

0 


Shuttle Antenna Gain 

- 4.0 dB 


Received Power 

-165.6 dBW 

At Antenna Terminals 

System Noise Temperature 
Boltzmann's Constant 

28.5 dBK 

-228.6 dB-W/K/Hz 

T = 705'-'’K (at antenna 

^ terminals) 

Noise Spectral Density 

-200.1 dB-W/Hz 


C/Nq 

34.5 dB-Hz 


Required C/Nq (Theoretical) 

28.5 dB-Hz 

= ' 5 ” 

Implementation Loss 

1.5 dB 


Required C/Nq 

29.5 dB-Hz 


Link Margin 

5.0 dB 



r 
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Table 6. Baseline Link Budget for Carrier Tracking, LI -C/A 


..ink LI -C/A, Carrier 
Upper Antenna 

Orbital Geometry: Path Tangential to Orbit 


P arameter Value Explanation 


Satellite EIRP 

26.8 

Space Loss 

-184.6 

Pointing Loss 

- 0.4 

Polarization Loss 

- 0.4 

Atmospheric Loss 

0 

Shuttle Antenna Gain 

- 4.0 

Received Power 

-162.6 

System Noise Temperature 

28.5 

Boltzmann'? “onstant 

-228.6 

:ioise Spectral Density 

-200.1 

C/Nq 

37.5 

Required C/N^ (Theoretical) 

28.0 

Implementatior Loss 

1.5 

Required C/Nq 

29.5 

Link Margin 

8.0 


dBW EOE. Spec, CID-SV-lOH 

dB Path Tangent to Orbit 

dB 
dB 

dB 

dBW At Antenna Terminals 

antenna 

terminals) 

dP-W/K/Hz 

dB-W/Hz 

dB-Hz 

^Mitter - ’5'’ 
dB 

dB-Hz 

dB 
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Table 7. Baseline Link Budget for Carrier Tracking, L2-P 


Link L2-P, Carrier 
Upper Antenna 

Orbical Geometry: Path Tangential to Orbit 


Parameter 

Value 

Explanation 

Satellite EIRP 

19.1 dBW 

EOE, Spec. CID-SV-lOH 

Space '.OSS 

-182.4 dB 

Path Tangent to Orbit 

Pointing Loss 

- 0.4 dB 


Polarization Loss 

- 0.4 dB 


Atmospheric Loss 

0 


Shuttle Antenna Gain 

- 4.0 dB 


Received Power 

-168.1 dBW 


System Noise Temperature 
Boltzmann's Constant 

28.5 dBK 
-228.6 dB-W/K/Hz 

T = 705®K (at antenna 

^ terminals) 

Noise Spectral Density 

-200.1 dB-W/Hz 


C/Ng 

32.0 dB-Hz 


Reqjired C/Nq (Theoretical) 

28.0 dB-Hz 

'’Jitter " 

Implementation Loss 

1.5 dB 


Required C/Nq 

29.5 dB-Hz 


Link Margin 

2.5 dB 
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Table 8. Baseline Link Budget for Carrier Tracking, L2-C/A 


Link 12-C/A, Carrier 
Upper Antenna 

Orbital Geometry: Path Tangential to Orbit 


Parameter 

Value 

Explanation 

Satellite EIRP 

19.1 

dBW 

EOE, Spec. CID-SV-lOH 

Space Loss 

-182.4 

dB 

Path Tangent to Orbit 

Pointing Loss 

- 0.4 

dB 


Polarization Loss 

- 0.4 

dB 


Atmospheric Loss 

0 



Shuttle Antenna Gain 

- 4.0 

dB 


Received Power 

-168.1 

dBW 


System Noise Temperature 
Boltzmann's Constant 

28.5 

-228.6 

dBK 

dB/W/K/Hz 

T = 705°K (at antenna 

^ terminals) 

Noise Spectral Density 

-200.1 

dB-W/Hz 


C/No 

32.0 

dB-Hz 


Required C/N^ (Theoretical) 

28.0 

dB-Hz 

^\litter ^ 

Implementation Loss 

1.5 

dB 


Required C/N^ 

29.5 

dB-Kz 


Link Margin 

2.5 

dB 
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Table 9. Baseline Link Budget for Data Detection, Ll-P 

Link Ll-P, Data Detection 
Upper Antenna 

Orbital Geometry: Path Tangential to Orbit 


Parameter 

Value 

Explanation 

Satellite EIRP 

23.8 dBW 

EOE, Spec. CID-SV-lOH 

Space Loss 

-184.5 dB 

Path Tangent to Orbit 

Pointing Loss 

- 0.4 dB 


Polarization Loss 

- 0.4 dB 


Atmospheric Loss 

0 


Shuttle Antenna Gain 

- 4.0 dB 


Received Power 

-165.6 dBW 

At Antenna Terminals 

System Noise Temperature 

28.5 dBK 

T = 705''’K (at antenna 

^ terminals) 

Boltzmann's Constant 

-228.6 dB-W/K/Hz 

Noise Spectral Density 

-200.1 dB-W/Hz 


C/Nq 

34.5 dB-Hz 


Required C/Nq (Theoretical) 

26.6 dB-Hz 

BER = 10'^ 

Implementation Loss 

2.0 dB 


Required C/Nq 

28.6 dB 


Link Margin 

5.9 dB 
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Table 10. Baseline Link Budget for Data Detection, LI -C/A 


Link LI -C/A, Data Detection 
Upper Antenna 

Orbital Geometry: Path Tangential to Orbit 


Parameter 

Value 

Explanation 

Satellite EIRP 

26.8 

dBW 

EOE, Spec. CID-SV-lOH 

Space Loss 

-184.6 

dB 


Pointing Loss 

- 0.4 

dB 


Polarization Loss 

- 0.4 

dB 


Atmospheric Loss 

0 



Shuttle Antenna Gain 

- 4.0 

dB 


Received Power 

-162.6 

dBW 

At Antenna Terminals 

System Noise Temperature 

28.5 

dBK 

T = 705®K (at antenna 

^ terminals) 

Boltzmann's Constant 

-228.6 

dB-W/K/Hz 

Noise Spectral Density 

-200.1 

dB-W/Hz 


C/Nq 

37.5 

dB-Hz 


Required C/Nq (Theoretical) 

26.6 

dB-Hz 

BER = 10"^ 

Implementation Loss 

2.0 

dB 


Required C/Nq 

28.6 

dB-Hz 


Link Margin 

8.9 

dB 



f 





14 


Table 11. Baseline Link Budget for Data Detection, L2-P 


Link L?-P, Data Detection 
Upper Antenna 

Orbital Geometry: Path Tangential to Orbit 


Parameter 

Value 

Explanation 

Satellite EIRP 

19.1 dBW 

EOE, Spec. CID-SV-lOH 

Space Loss 

-182.4 dB 

Path Tangent to Orbit 

Pointing Loss 

- 0.4 dB 


Polarization Loss 

- 0.4 dB 


Atmospheric Loss 

0 


Shuttle Antenna Gain 

- 4.0 dB 


Received Power 

-168.1 dBW 

At Antenna Terminals 

System Noise Temperature 
Boltzmann's Constant 

28.5 dBK 
-228.6 dB-W/K/Hz 

T, = 705*K (at antenna 
terminals) 

Noise Spectral Density 

-200.1 dB-W/Hz 


C/No 

32.0 dB-Hz 


Required C/Nq (Theoretical) 

26.6 dB-Hz 

BER = 10"^ 

Implementation Loss 

2.0 dB 


Required C/Nq 

28.6 dB-Hz 


Link Margin 

3.4 dB 




Table 12. Baselink Link Budget for Data Detection, L2>C/A 


Link L2-C/A, Data Detection 
Upper Antenna 

Orbital Geometry: Path Tangential to Orbit 


Parameter 

Value 

Explanation 

Satellite EIRP 

19.1 

dBW 

EOE, Spec. CID-SV-lOH 

Space Loss 

-182.4 

dB 

Path Tangent to Orbit 

Pointing Loss 

- 0.4 

dB 


Polarization Loss 

- 0.4 

dB 


Atmospheric Loss 

0 



Shuttle Antenna Gain 

- 4.0 

dB 


Received Power 

-168.1 

dBW 

At Antenna Terminals 

System Noise Temperature 

28.5 

dBK 

* 705'^K (at antenna 




terminals) 

Boltzmann's Constant 

-228.6 

dB-W/K/Hz 


Noise Spectral Density 

-200.1 

dB-W/Hz 


C/Nq 

32.0 

dB-Hz 


Required C/Nq (Theoretical) 

26.6 

dB-Hz 

BER = 10"^ 

Implementation Loss 

2.0 

dB 


Required C/Nq 

28.6 

dB-Hz 



Link Margin 


3.4 dB 



Table 13. Baseline System Configuration, On-Orbit 
Link Margin Summary 



Link Margin by Function 


Link 

Ranging 

m) 

Carrier 

(dBl 

Data 

(BB) 

Lower Antenna* 

LI P 

3.3 

5.4 

6.3 

C/A 

6.3 

8.4 

9.3 

L2 P 

2.9 

2.9 

3.8 

C/A 

2.9 

2.9 

3.8 

Upper Antenna 

LI P 

2.9 

5.0 

5.9 

C/A 

5.9 

8.0 

8.9 

L2 P 

2.5 

2.5 

3.4 

C/A 

2.5 

2.5 

3.4 


*0V102 lower antenna margins are 1.3 dB lower. 





This geometry, called LOS tangent-to-orbit, Is considered a 
nominal baseline case. Some links, such as the overhead viewing links, 
will have less space loss: and other links, such as those where the 
Orbiter may be in higher orbits and viewing tangent to the earth (or 
earth's atmosphere), will have greater loss. 

Pointing Loss 

The pointing loss value of -0.4 dB occurs because the GPS satellite 
antenna coverage is designed to provide maximum radiated power at the ed^e 
of the e«»rth. Since the Orbiter is some 200 miles above the edge of tha 
earth, there is a reduction in EIRP. The GPS satellite antenna pattern 
for FSV #1 has a gain slope of approximately 0.5 dB/degree. From Figure 2, 
the angle seen by the GPS satellite antenna is 0.72 degree for the 200-mile 
Shuttle orbit. Thus, the EIRP loss is equal to 0.5 JB/degree x j.72 degree 
* 0.4 dB. 

Polarization Loss 

The polarization loss occurs because the GP'!‘ antenna and the 
Orbiter GPS antenna are not perfectly circular-polarized. GPS satellite 
antenna test data indicates an axial ratio of approximately 1.2 dB at the 
edge of the earth. Assuming a Shuttle antenna axial ratio of 4 dB and 
referring to Figure 3, which gives polarization loss between two ellip- 
tically polarized antennas, the loss is seen to be approximately 0.4 dB. 

Antenna Gain ^ 

The -4 dB gain is with respect to perfect RHCP and has been choseii 
by Rockwell to give the desired angular coverage. It is based on breadboard 
antenna measurements and is subject to change. It should be noted that a 
decrease in overall actual antenna gain can translate into a decrease in 
antenna coverage volume rather than a decrease in link margin. 

Circuit Loss 

The circuit losses used for the link budget calculations have been 
provided by Rockwell and are given in Table 14. Two types of losses are 
given: the loss from the antenna tc the preamplifier and the loss from 

the preamplif'*er to the receiver. The latter loss includes a 6 dB power 
splitting loss for adding a third-string receiver. Even though the base- 
line design is a two-string system, it is cost effective to install and 










20 


origin 
OF POOR 


pagk is 
qUAU'I^ 



Figure 3. Polarization Loss Between Two Elliptically Polarized Antennas 
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connect the necessary power splitters for a three-string system. 


^"-'-'--....Circuit 

Antenna Circuit (L,) 
(dB) ^ 

PA/RPA Circuit (Lo) 

Antenna^'-"'--....^ ^ 

(dB) 

Upper Antenna 

1.66 

11.22 

Lower Antenna 

1.36 

11.88 

(OV 102) 

(2.42) 

(11.88) 


Table 14. RF Circuit Loss 


System Noise Temperature 

The signal -to-noise ratio required for all GPS receiver processing 
functions is proportional to C/kT^ (C/Nq) where C is the available signal 
power, k is Boltzmann's constant, and T^ is the effective noise temperature. 
Since the Shuttle GPS receiving system is comprised of a cascade of ele- 
ments, it is necessary to derive the effective noise temperature T^ for 
the desired point for definition of C/KT^. 

The Shuttle GPS system is modeled in the block diagram shown in 
Figure 4. As will be shown, it is inconsequential as to which point is 
picked as a reference point for definition of C/Nq since the ratio is con- 
stant, i.e., independent of the point in the circuit. A generalized model 
of a cascaded system is used to illustrate this point. This model is 
shown in Figure 5. The effective noise temperature at the input to the 
cascade. T^, is given by 

To T^ 

‘e 'l G^ G^G^ 

and the input signal power 
temperature at point (J) , 
to the first box will have 
Thus, 


+ - - + 
G j G2G^ 


is C. Now, to find the system ettective noise 
T ' , we note that noise at from the input 
been increased by the gain of this box, . 


G.T 

^ 1 T 

^i^M ^ ^ Go ^ Go?,; ' 
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Figure 5. General Mode of System Noise Temperature for 

Cascaded Elements 
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Similarly, the signal power at (|) is increased by 6^ so that 


C » 6^C 


and the signal-to-noise ratio is given by 

C 


G^C 


KT, 





which, of course, is the same as the signal-to-noise ratio at the input 
to the cascade. The same reasoning can be applied to show that the 
signal-to-noise ratio at (3) is also the same by multiplying both signal 
and noise at (i) by gain G^. The noise temperature at (s) is given by 


T,“ = TjGjGj * T^Gj . Tj * 


The generalized model of Figure 5 is now applied to the Shuttle 
GPS system. Box 1 is taken as the antenna with temperature T^, box 2 is 
the loss between the antenna and preamp, L-| , box 3 is the preamp, box 4 
is the loss between the preamp and the receiver, L 2 , and box 5 is the 
receiver. The equivalent noise temperature of an active element of noise 
figure NF and gain 6 is given by 

T = (NF-1) Tg 

where Tg is the reference temperature (290°K). The noise temperature of 
a passive element of loss L and gain G = 1/L is given by 

T = (L-l)Tg 


Thus the equivalent noise temperature of the GPS system at the antenna 
terminal (Tg') is given by 




(NFpj-l)To (Lj-DTo ^ (NfRPA-')T( 


1/L 


1 '/L, GpA ’/h 


l/Lj 


The system noise temperature and the signal power in the link 
budgets which follow are referenced to the GPS antenna terminals. There- 
fore, the system signal power- to-noise density ratio is written as 


t 
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; 

c 


kT, 


T^Mh-l)^L,{NFp^-l)^L^ 


H4(NFppA-l) 

®PA ®PA 


The system equivalent noise temperatures for the GPS system have 
been calculated for the various antenna and loss combinations using the 
parameters 

NFp^ = 3.0 dB 
GpA * 37 dB 
NF,„ - 13 dB 

,l 2 “ values given in Table 14. 

Table 15 presents a summary of the system equivalent noise temper- 
atures, both in degrees kelvin and in dB-®K. 


Table 15. Summary of Effective System Noise Temperature 
(at antenna terminals) 


Antenna 

Effective Temperature (includes antenna 
temperature = 125°K) 


Degrees kelvin 

dB-®K 

Upper 

705 

28.5 

Lower 

650 

28.1 

(0V102) 

876 

29.4 


It should be noted that 0V102 lower antenna circuit noise is higher than 
the other two cases shown in Table 15 due to the diplexer loss as a result 
of sharing the antenna feed cable with the S-band antenna. 

Required C/Nq and Implementation Loss 

The required C/Nq for ranging is based on a range measurement 
accuracy of a = 5 feet for Ll-P and a = 7 feet for L2-P. The L2 



requirement is relaxed because L2 Is needed on1> for Ionospheric correc- 
tion and its measurement is effectively averaged over several samples 
in the correction process. The derivation of e 29.6/27.5 dB, as well 
as the 2 dB implementation 1oss» is documented in [1]. In addition, 
the carrier tracking requirement of 28 dB is based on an allowable RMS 

phase jitter of 15**, and the data detection requirement of 26.6 dB is 
-5 

based on a 10 BER for the 50 bps data. These parameters are also 
documented in [1]. 
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3.0 SHUTTLE GPS PREAMPLIFIFR DESIGN 


The Shuttle GPS system requires an RF preamplifier mounted close 
to each antenna to control the overall system noise figure. This is due 
to the fact that there is a long length of feedline between each antenna 
and the GPS receiver processor assembly (RPA), in addition to other loss 
elements such as RF switches. An equivalent block diagram of this system 
is given in Figure 6. The equivalent noise temperature of the system, 
referenced to the input of the RPA is given by 



(L,-l)To* 


(NFpj^-DTg (4-DTo 

!7t; (i/L,)Gp^ 


(NF^3 

-nTqt^- 


Thus, it is obvious that in addition to minimizing the antenna feed loss 
it is most important to minimize the preamp noise figure, NFp^, while at 
the same time maximizing the gain so as to minimize the noise contributions 
of loss L 2 and the RPA front end. 



Figure 6. Model of GPS Receiving System for Calculation of 

Equivalent Noise Temperature 


A recommended candidate Shuttle GPS preamplifier design is shown 
♦ in block diagram form in Figure 7. Each of the elements in Figure 7 will 

be discussed for the purpose of evaluating the feasible performance for 
this GPS preamplifier. An alternate preamplifier configuration is shown 
in Figure 8. This configuration functionally performs the same; however, 
t it does not offer the redundancy advantage that the configuration of 

Figure 7 offers. The advantage of figure 8 is that it eliminates one 
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Amplifier DC On/Off Must Be 
Coordinated with Switch Position 


Figure 7. Recommended Redundant Preamplifier Configuration 



Ampl ifier 


Figure 8. Recommended Nonredundant Preamplifier Configuration 
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amplifier, the switch, and the summer. The major disadvantage of the 
Figure 8 preamplifier configuration is the lack of redundancy. If the 
amplifier fails, the antenna to which the preamplifier is connected is 
no longer useable for the Shuttle GPS navigation system. Of course, in 
the case of the redundant preamplifier, if the antenna itself fails, the 
preamplifier redundancy is of no use in providing the signal reception. 

One argument, however, in favor of the redundant preamplifier configuration 
is that the MTBF of the antenna is several orders of magnitude greater 
than an active electronic circuit. This is the reason why communication 
satellites, and most scientific spacecraft, almost always use redundant 
transponders (and preamplifiers, when preamplifiers are used), even when 
there is a single antenna. 

Yet another preamplifier configuration is shown in Figure 9. This 
is the configuration used in the Phase 1 GPS developmental receivers. It 
offers the advantage of dual redundancy without having to utilize a switch. 
The redundancy gained is not as powerful as the redundancy provided by the 
Figure 7 configuration, since in the case of Figure 9 one preamplifier 
channel is dedicated to the LI frequency and the other channel is dedicated 
to the L2 frequency. Thus, a failure of one channel would constrain opera- 
tion (with the associated antenna) to the unaffected frequency. A further 
drawback to this configuration is that the diplexer used to separate the 
L1/L2 channels adds significantly to the cost and size of the preamplifier. 

The broadband amplifier shown in Figure 7 is the one element which 
most strongly determines the performance of the preamplifier. It is desir- 
able for this amplifier to have as low a noise figure as possible, along 
with high gain, but consistent with reasonable cost. In order to determine 
reasonable performance specification requirements for the Shuttle GPS pre- 
amplifier, a survey has been taken of some of the off-the-shelf commercial 
low noise GaAs, FET, and bi -polar amplifiers. A summary of amplifier 
characteristics for amplifiers covering a wide frequency range is given 
in Table 16. The circled line represents the amplifier most closely satis- 
fying the Shuttle GPS requirements. The two parameters of greatest signi- 
ficance are the noise figure of 2.5 dB and the 38 dB gain. It must be 
stressed that these values are the minimum values over the frequency range 
of the amplifier. Typical variation with frequency is illustrated in 
Figure 10. The preamplifier performance is not entirely determined by 



Table 16. Typical 


FREQ. 

RANGE 

(MH{) 

NOISE 

FIGURE 

(dBI 

MAX. 

GAIN 

(dBI 

MIN. 

9Sai250 

2.2 

24 

loaaiooo 

22 

24 

1100-1400 

22 

30 

llOSO-IBSO 

26 

38 

1430-1540 

26 

24 

isoa2ooo 

3.0 

50 

17sai850 

2.8 

33 

1400-2400 

26 

M 

1400-2400 

20 

28 

1700-2400 

3.5 

28 

2100-2400 

3.3 

28 

2200-2300 

3.1 

28 

2200-2300 

3.3 

28 

2100-2400 

4.3 

28 

2400-2700 

25 

29 

2S0a2S80 

1.8 

30 

2500 2580 

2.2 

30 

2700-2900 

3.5 

24 

2700-2900 

5.5 

4i4 

2700-3100 

3.7 

24 

310a3500 

4.3 

30 

3100 3500 

6.0 

30 


f 
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Bl-Polar and GaAs FET Low Noise Amplifiers 


fo 

(dBrnl GAIN INTERCEPT 

MIN 9 FLAT POINT 

I de NESS (dBm) 

COMP (idB) TYP 

♦ 5 as *16 

+ S 0.? +16 

♦10 1.0 ^21 


♦10 

1.0 

♦21 

2.0 

2.0 

401 

^ 8 

23 

♦15 

2.0 

1.1 

23 

415 

1.0 

♦26 

2.0 

ZO 

65 

410 

0.3 

+21 

2.0 

2.0 

35 

410 

1.0 

♦21 

2.0 

2.0 

36 

413 

1.0 

♦24 

2.0 

10 

45 

410 

0.8 

♦21 

2.5 

1.1 

» 

4 9 

0.5 

♦20 

2.5 

1.5 

35 

4 9 

0.3 

♦20 

2.5 

1.3 

35 

* 9 

0.3 

♦20 

1.3 

1.3 

35 

4 9 

0.5 

+20 

1.5 

1.6 

35 

410 

0.8 

♦21 

2.0 

2.0 

40 

♦ 9 

0.3 

♦20 

1.3 

1.4 

80 

♦ 9 

0.3 

+20 

1.3 

1.4 

80 

411 

0.5 

+22 

2.0 

2.0 

40 

412 

0.5 

♦23 

2.0 

2.0 

45 

410 

0.5 

♦21 

2.0 

1.7 

40 

410 

0.5 

♦21 

2.0 

2.0 

50 

410 

0.5 

♦21 

2.0 

2.0 

50 


VSWR 

MAX OC POWER 
4ISV0C9 


JN 

OUT 

(mA)NOM 

2.0 

1.5 

23 

2.0 

1.5 

23 

2.0 

£0 

30 
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the low noise amplifier performance due to the additional elements. The 
effect of these elements Is now discussed. 

The purpose of the EMI filter is to prevent strong unwanted signals 
that are outside the GPS frequency bands from saturating the preamplifier. 
The simplest form of filter is the bandpass filter that passes both LI 
and L2 frequencies, as shown in Figure 11a. The disadvantage of this 
filter is that interfering signals that fall in the 300 MHz spectral 
region between LI and L2 can potentially saturate the preamplifier. To 
eliminate this problem, the filter would need to be a double tuned filter, 
peaked at the LI and L2 frequency of 1575 MHz and 1227 MHz, respectively, 
as shown in Figure lib. The insertion loss for this filter will be some- 
what greater than that for the bandpass filter case. 

Since the filter precedes the low noise amplifier, as does ihe 
switch, insertion loss from these sources directly cegr/'Hes the preampli- 
fier overall noise figure. Thus, a tradeoff exists between filter selec- 
tivity, which is obtained by adding more sections and thus increasing the 
insertion loss, and overall preamplifier noise figure. TK' selectivity 
of the bandpass filter as a function of the number of sections is illu- 
strated in Figure 12, and the insertion loss of the filter as a function 
of the number of sections is given by 


IL 


T I.... _ (Loss Constant) (Number of Sections + 1/2) « o 

Insertion Loss (.ertenf 'J-JOS * 


Thus, for a three-section filter, a reasonable compromise that has a 3 dB 
bandwidth of 400 MHz (centered on — + 1227= 1400 MHz), the inser- 
tion loss is found to be given by 

' feo/ 1400)1 00 
= 0.35 dB . 

From Figure 12, it is seen that the filter response is down 20 dB at fre- 
quencies of 680 MHz and 2120 MHz. Since this roll-off is probably slower 
than the amplifier roll-off, a four-section filter is chosen. This gives 
an insertion loss of 0.39 dB, and a response 20 dB down at approximately 
880 MHz and 1920 MHz. 



1227 MHz 


1575 MHz 


Figure 11 


. Bandpass Filter Shape with No EMI 
Between LI and L2 



b. Bandpass Filter Shape with EMI Threat 
Between LI and L2 


. GPS Preamplifier' Bandpass Filter Shapes 
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Figure 12. Bandpass Filter Attenuation Characteristics 
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The RF switch shown in F1<r v 7 is implemented as a magnetic 
latching switch. This means that t , switch is normally deenergized 
and in one of the two stable positions. To switch, a short pulse is 
applied to the appropriate switching coil. This pulse is typically 
approximately 15 milliseconds in duration and about 2 watts. The outline 
drawing for a typical low-loss coaxial switch is shown in Figure 13. The 
significant performance parameter for the switch is the insertion loss 
of 0.2 dB. 


rvf 



Figure 13. Outline Drawing for Typical Low Loss 
Coaxial Microwave Switch 


The signal sunmer shown following the broadband amplifiers in 
Figure 7 is a 3 dB hybrid power summer. Use of a summer for this element 
means that the unused amplifier must be turned off (DC power turned off) 
to avoid a 3 dB degradation in preamplifier noise figure. To avoid turning 
off the unused amplifier, another switch can be substituted for the 3 dB 
hybrid. 

The overall noise temperature of the preamplifier, considering the 
elements discussed above, is given by 

(NF-l)T. (U-1)T» 

TpA = (L^ - 1 )Tq + “T7T^ * 

where L^ = combined loss of input selection switch and input filter 
NF = noise figure of low noise amplifier 
G * gain of low noise amplifier 
L 2 = loss of output switch or 3 dB hybrid . 



The effective noise figure of the preamplifier is given by 




1 


NFL^ + 


(L2-1)^ 

G 


From the survey of available commercial components, the following param- 
eter values are appropriate: 

NF = 2.4 dB 
G = 38 dB 

= 0.2+ 0.4 dB = 0.6 dB 

L 2 = 0.2 dB (switch) = 3.2 dB (hybrid) . 

Thus, the overall noise figure for the preamplifier is found to be 3 dB. 

It should be noted that the choice of either the switch or the hybrid for 
the output selection affects the overall noise figure so insignificantly 
that it is not a factor in the preamplifier design configuration decision 
nor in the overall system performance. 

Since the preamplifier noise figure is so critical to system per- 
formance, techniques for testing this parameter are discussed here. 

Using a test set-up such as shown in Figure 14, the receiver is 
adjusted to the desired NF measurement frequency. The noise figure meter 
input selector is set to the receiver IF output frequency. 

The variable attenuator is used to optimize the input level to 
the NF meter. 

Precautions should be taken when using a gas tube noise source to 
assure that large voltage "spikes" will not be applied to the amplifier 
input, with resultant damage or erroneous results. Also, where necessary, 
the noise figure reading obtained should be corrected for noise source 
variances (from nominal specified excess noise ratio) and for second stage 
(receiver or post amplifier) contributions. When such precautions and 
corrections are adhered to, good correlation is obtained between various 
commercial noise figure measuring equipments. 
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Since the method of selecting between redundant preamplifiers is 
an electromechanical switch, as shown in Figure 7, it is of interest to 
consider alternate methods of selecting the preamplifier. One such alter- 
nate, which at first appears attractive, is a 3 dB hybrid power splitter. 

The analysis which follows shows that this can lead to a 3 dB signal-to- 
noise degradation and is thus not a satisfactory choice. Furthermore, 
an analysis is presented which shows that a switch rather than a hybrid 
should be used at the receiver to select between the upper and lower 
antennas. 

The simplest way to provide dual redundancy for the preamplifiers 
is to parallel them by using 3 dB hybrid power splitters/summers, as shown 
in Figure 15. The equivalent noise model is shown in Figure 16. The ampli- 
tude of the output (input to ideal summer) of the upper arm is given by 

/2 vT" 

and the output of the lower arm is 


SO that the output of the summer is given by 


VgCt) = g 


ni(t)+ n,(t) 
(tl + -! — 


s(t) 


2 2 2 

The signal power is g s(t) =g P, and the noise power is given by 


N = 


^2 


g [H^(t) + 


= g 


n,(t)^ n„(t)2 

-J + -£ 


;2 


. 2 

-go. 
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Figure 15. Shuttle/GPS RF System Block Diagram 
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2 

SO that the output signal-to-noise ratio is S/N*P/o , which is the same 
signal-to-noise ratio found at the output of the single preamplifier without 
the power splitter. Thus, there is no loss in SNR in parallelling two 
preamplifiers with hybrid power splitters/summers. However, this is not 
the case if one of the preamplifiers fails. In this case, the output of 
the upper arm is still 

-s(t) . 


but the output of the lower arm is zero, by setting the gain g = 0. Thus, 
the output of the summer is 


Vo(t) 


gs(t) ^ 


and the signal power is 

s(t)^ , 

4 4 


and the noise power is given by 


N 


g^n^(t)^ 

2 



» 


2 

so that the output signal-to-noise ratio is S/N = P/2o . Thus, if one of 
the parallel preamplifiers fails, there is a 3 dB loss in SNR due to the 
3 dB signal loss in the input power splitter. This loss is avoided by 
replacing that splitter with a SPOT switch. 

The signal-to-noise ratio that results from feeding one receiver 
from two antenna/preamplifiers via a hybrid summer can be found by 
rearranging the model in Figure 16 into the model shown in Figure 17. 
There is no signal in the lower arm since only one antenna at a time 
picks up the desired GPS signal. The output of the upper arm is 


/r >fT 


» 


and the output of the lower arm is 





I 




t 


I 






43 


gn9(t) 


vT 

so that the sunmer output Is 

fl^+ g 


n^(t) n2(tri 




Jl 


The output signal -to -noise ratio is 


SNR 


n^(t) ^ n 2 (tr[‘ 


2o 


2 • 


Thus, there is a 3 dB loss in signal-to-noise ratio if a hybrid summer 
is used to sum the upper and lower antenna inputs at the RPA. For this 
reason, a switch must be used for the antenna management function. 


♦ 


♦ 
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4.0 SHUTTLE GPS RECEIVER ARCHITECTURE 


There are two basic GPS receiver architectures that have been 
considered for the Shuttle GPS receiver. They are the parallel channel 
receiver and the sequential receiver. Examples of the parallel channel 
receiver are the X set and the Collin's GDM set. A functional block dia- 
gram of the X set is shown in Figure 18. Examples of sequential receivers 
are the manpack receiver and the GPSPAC receiver. A functional block dia- 
gram of the GPSPAC is shown in Figure 19. 

The major generic difference between the parallel channel receiver 
and the sequential receiver is a simple one. The sequential receiver 
time multiplexes one channel among four or more satellite signals in a 
serial manner. In the parallel channel receiver, four or more parallel 
hardware channels simultaneously process the satellite signals. Somewhat 
better performance and platform dynamics capability are the advantages 
gained by the parallel channel set. However, this is at the expense of 
almost four times the hardware complexity of the sequential receiver. 

This translates into greater volume, weight, power consumption and cost 
for the parallel channel set. Of great significance in the decision pro- 
cess for choosing the baseline Shuttle GPS receiver is that no existing 
parallel receiver would fit within the volume constraints of the avionics 
equipment bay. 

On the other hand, the size of the existing sequential receivers 
was found to be compatible with the Shuttle volume constraints. Further- 
more, analysis has shown that the navigation performance of the sequential 
receiver meets the requirements of a Shuttle GPS navigation system. Anal- 
ysis has also shown that the sequential receiver will track the Shuttle 
signal dynamics. Consequently, a sequential receiver architecture has 
been chosen for the baseline Shuttle GPS receiver. A second channel for 
the sequential receiver is recorranended to provide faster acquisition (time 
to first fix) and a measure of redundancy. 

It should be noted that this choice should not rule out considera- 
tion of future parallel channel receivers that take advantage of technology 
to meet the strict volume constraints of the Shuttle. 

The remainder of this section presents a discussion of some of the 
more detailed trade-offs that can be made within the sequential receiver 
architecture; first, however, a generalized review of the sequential 
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Figure 18. Functional Block Diagram of X Set Parallel Channel GPS Receiver 
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Figure 19. 6PSPAC Receiver Functional Block Diagram 












receiver concept is presented. 

A simplified functional block diagram of a single-channel 
sequential receiver is shown in Figure 20. The salient feature of this 
type of receiver is a single RF/IF chain that is time-multiplexed among 
satellite measurements. Time-multiplexing is also utilized for tracking 
loops, and the other baseband functions such as data detection, acquisi- 
tion and synchronization detection, etc. A typical sequence of operations 
for this type of receiver is shown in Figure 21. The basic time element 
of this approach is the frame which is divided into satellite slots. 

Each satellite slot is divided into an LI and an L2 slot and these slots 
are divided into the receiver process intervals. Typically, two seconds 
are sufficient for an LI ot L2 slot so that, for processing four satellite 
signals, a 16-second frame is required. If two parallel receiver channels 
are used so that an LI and L2 can be processed simultaneously, the frame 
time can be halved, i.e., an 8-second frame. However, it has been shown 
that a shorter frame time does not provide any beneficial navigation 
performance for the Shuttle GPS system. A parallel L2 channel would 
increase receiver complexity and cost and is thus a poor trade-off. It 
should be noted, however, that parallel channels in the sequential 
receiver may be desirable from the redundancy and initial acquisition 
time (time to first fix) viewpoint. 

The first receiver-processing interval is used to search and 
acqu're the P code. This involves prepositioning, or slewing, the coder 
to the predicted code phase for the particular satellite to be acquired. 
This phase information is calculated and stored by the receiver processor 
(microcomputer) for each satellite on a frame -to-frame basis. Once the 
coder is slewed to the proper phase, it goes into a limited automatic 
search mode wherein it searches out the few code chips of uncertainty 
(see Section 6.0 for code acquisition discussion). At this same time, 
the receiver AFC is enabled. As the correct code phase is reached, the 
code search stops and the signal -to-noise ratio increases due to correla- 
tion so that the AFC loop starts to pull the receiver VCO to the correct 
frequency. 

In the second receiver-processing slot, the code loop is enabled 
and starts to pull into phase lock, thereby decreasing the correlation 
error to a small fraction of a chip. The AFC pulls in the frequency more 
rapidly due to the increase in S/N. 
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Figure 21. Sequential Receiver Timing 
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During the third receiver-processing slot, the Costas loop is 
enabled and starts to pull in. In addition, the code loop is disabled 
since, by this time, the loop has pulled the code phase error into an 
acceptable error. Since the PN code is coherent with the carrier, the 
Costas loop VCO can be scaled by the frequency ratio (1575.42/10.23) and 
used as a clock source to drive the coder. 

Finally, in the fourth slot, with the Costas loop tracking in 
steady state, the pseudo range and delta range (range rate from carrier 
times interval) measurements are made. 

The sequential receiver process is highly dependent on micropro- 
cessor control. Figure 22 illustrates the interfaces between the 
receiver microprocessor and the other receiver functions. In this par- 
ticular architectural scheme, the carrier-tracking function is implicit 
in the synthesizer wherein the phase and frequency of an appropriate 
local oscillator (LO) is adjusted by the rate multiplier/incremental 
phase modulator (RM/IPM). The RM/IPM is, in turn, controlled by the 
Costas loop error signal developed in the microprocessor. All the 
tracking error signals are developed in the microprocessor by appropriate 
processing of the I and Q signals. The I and Q signals are at baseband 
and the A-to-D converter that converts these signals to digital words is 
the major interface between the analog and digital portion of the receiver. 
It should be pointed out that it may not be practical to implement the 
RM/IPM function digitally for the Shuttle GPS receiver due to the high 
doppler offsets. An analog VCO might be needed. Whether a digital VCO 
or analog VCO is used, the loop filtering and error control are via the 
microprocessor. 

The RF/IF frequency scheme has particularly strong impact on over- 
all receiver architecture due to the impact on synthesizer complexity. 

The general issue is whether single conversion, double conversion or even 
triple conversion from RF to IF to baseband should be used. From an over- 
all receiver complexity consideration, single conversion provides the 
least complex receiver. On the other hand, this type of receiver presents 
the greatest design risk due to the large gain required at a single fre- 
quency. In the discussions of frequency conversion trade-offs that fol- 
low, translation to baseband (I and Q channels) via coherent LO is not 
treated as a frequency conversion in the traditional sense but, rather, 




Figure 22. Digital GPS Receiver Interfaces 
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as part of the traditional demodulation process. 

The frequency conversion plan for a single-conversion GPS receiver 
that is designed to minimize volume and cost is shown in Figure 23. Both 
the first LO and the coherent reference for translation to I and Q are 
derived from the same VCO which is part of the carrier tracking loop. 

This receiver architecture is simple and straightforward and has a simi- 
larly simple frequency synthesizer; however, this receiver requires 
approximately 120 dB to 140 dB gain at the single IF frequency of 
71.61 MHz. This much gain at a single IF frequency is counter to tra- 
ditional receiver design practices. Careful evaluation of the detailed 
receiver design is required to assess the design risk. 

At the opposite extreme of receiver design is the triple-conversion 
approach to GPS receiver design, as shown in Figure 24. This frequency 
plan distributes the gain at three IF frequencies; it uses three local 
oscillators in addition to the coherent translation to I and Q baseband 
channels. This receiver architecture has been used in the manpack and 
GPSPAC receivers with good performance results. As might be expected, 
however, it results in a complex frequency synthesizer. 

A potential compromise between the single- and triple-conversion 
GPS receiver architectures is the double-conversion receiver, as shown 
in Figure 25. The receiver design is of moderate complexity and offers 
the advantage of distributing receiver gain at two IF frequencies. A 
possible frequency synthesizer plan for this receiver is shown in 
Figure 26. The plans in Figures 25 and 26 are presented as possibilities. 
Detailed studies of potential spurs and image frequency problems are 
necessary before an actual design of the double-conversion receiver 
can be reconmiended. 

As a result of the investigation of the basic trade-offs amor.) 
single-, double-, and triple-conversion receiver architectures, it appears 
that double-conversion represents a good choice. This, of course, must 
be weighed against any existing receiver designs which have proved suc- 
cessful and which are presented for Shuttle application on a competitive 
performance, volume, and cost basis. 


PN Code 
(Baseband) 






Detector 







Kigure 25. Double Conversion GPS Receiver Frequency Conversion Plan 










5.0 SHUTTLE GPS EMI/RFI ISSUES 


This section discusses the analyses of Shuttle emitter RF 
harmonics and intermods interfering with GPS. It also examines the 
antenna isolation required between the TACAN antenna and GPS antenna 
to minimize degradation to the GPS navigation system. The conclusions 
of both analyses are that there is no interference problem. 

5.1 Evaluation of Intermodulation Interference 

Because the Shuttle has several RF emitters on board capable of 
radiating simultaneously, there is a possibility of intermods falling 
within the GPS frequency bands. Furthermore, there is also the possibil- 
ity of harmonics of the lower frequency emitters either falling on GPS 
or mixing with other emitters and falling on GPS. Since this problem is 
not amenable to simple analytical techniques due to the large number of 
possible frequency combinations, a JSC computer program is used. This 
program, operated by Art Reubens, is shown in a functional diagram in 
Figure 27. As can be seen, the basic inputs are transmitter characteris- 
tics, receiver characteristics and antenna characterisics. The program 
calculates all harmonics up to the fifth harmonic and all intermods up 
to the fifth order. After calculating the antenna coupling factors 
involved, the program searcnes the intermods and harmonics for ones that 
fall within the specified bandwidth, in this case, 30 MHz centered about 
Ll and L2 GPS frequencies. Furthermore, it only prints those intermods 
that fall above a specified minimum power level in this case of -130 dBm. 
A sample output sunwiary of the program is given in Table 17. 

A majority of the interference found can occur only at the launch 
pad since the intermods are a result of mixing Shuttle emissions with 
KSC C-Band and/or weather radars. Axiomatix recommends that this form 
of interference not be taken into account in the Shuttle/GPS design 
trade-offs. This is because, should this interfet nee actually exist, 
operational procedures can be used which will enable a launch pad GPS 
navigation fix to be obtained prior to launch. Furthermore, this inter- 
ference is relatively insensitive to GPS antenna locations. For obvious 
reasons, those interferences which are "not operationally possible" have 
not been considered further. The only interferences of possible concern 
are those labeled as "significant products." Each of these products. 
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Figure 27. Computer Evaluation of Intermodulation Interference Functional Diagram 
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Table n. Computer Evaluation of Intermodulation Interference, 

Sample Output 


Source 

Frequency 

Harmonic Beat Frequency 

Power Level 

TACAN Channel 35X 
Return to Landing 
Site at Kennedy 

1,059 MHz 

2nd 


UHF 

296.8 MHz 

3rd 




1,227.6 MHz 

-144.5 dBm 
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a total of 134, has been examined by Axiomatlx. All but four of them 
are a result of mixing at least two Orbiter TACAN emisslohs with anothe.' 
Orbiter emission, with each TACAN operating on a separate, specific 
channel. Thus Axiomatlx recotr^mends against considering these Interfer- 
ence sources In the Shuttle/GPS design trade-offs since It should be 
relatively easy to plan an operational “work-around" for obtaining a GPS 
navigation fix when two or more TACANs happen to be transmitting simul- 
taneously. The four cases which can be considered potential operational 
sources of Interference are summarized In Table 18. It should he noted 
that those four cases are all with the same basic set of emitters, the 
variables being antenna locations. Thus the frequency of the Intermods 
summarized In Table 18 Is constant, while thv amplitude varies as a func- 
tion of antenna geometry. The worst case Is number 2, with the Inter- 
fering Intemod having a power level of -112 dBm. The effect of this 
Intermod on the GPS signal processing Is to raise the effective thermal 
noise level. Since the 10 megachip PN correlation process spr» .»s the 
Intermod power Into an effective noise bandwidth of 10 MHz, the noise 
density of the "spread Intermod" Is then -112.9 - 10 log (10 MHz) = 

-182.9 dBm/Hz. The natural thermal noise density of nominally -171 dBm/ 
Hz so that the Interfering noise is 11.8 dB lower and effectively adds 
about 0.3 dB to the overall system noise. The LI and L2 P code links 
can easily accommodate this degradation. Furthermore, since the Intermod 
analysis Is a worst-case analysis, the actual Interference will probably 
not be as strong. Since the Intermod spectnan Is wider (approximately 
5 MHz) than the 1 megachip C/A code signal, the power spectrum density 
for the C/A Interference case Is -112.9 - 10 log (5 MHz) * -179.9 dBm/Hz 
which results in approximately 0.6 dB degradation. The C/A link can sup- 
port this degradation but, more Importantly, since the radar altimeter is 
used below 2500 feet altitude, the GPS receiver will be operating with 
the P links in order to have navigation accuracies compatible with the 
low altitudes. 

5.2 GPS/TACAN Antenna Isolation Requirements 
5.2.1 Summary 

The required isolation between the GPS antenna and the Shuttle 
TACAN antenna has been analyzed and found to be 28 dB. The analysis is 
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based on extrapolating the TACAN transmitter output spectrum measurements 
made by JSC. 


5.2.2 Discussion 


The requirements for GPS antenna isolation from the Orbiter TACAN 
antennas has been analyzed in light if the information contained in JSC 
report EE6-12/77-189, “TACAN Transmitter Pulse Characteristics." This 
antenna isolation is based on the desire to have the TACAN wideband modu- 
lation sideband energy falling in the GPS RF bandwidth to contribute less 
than 0.5 dB to the receiver overall noise figure. Since the TACAN PRF is 
relatively low (150 pps max) and the GPS precorrelator RF bandwidth is 
relatively wide (10 to 20 MHz) and separated from the TACAN center fre- 
quency by at least 77 MHz, the TACAN energy will appear as uncorrelated, 
wideband noise. 

The natural noise spectral density for the Orbiter GPS receiver 
is approximately -171 dBm/Hz. Thus, if the spectral density of the 
TACAN interference in the GPS bandwidth is limited to -181 dBm/Hz, the 
degradation to the total signal-to-noise ratio will oe less than 0.5 dB. 
Consequently, the difference between the actual TACAN interference spec- 
tral density and -181 dBm/Hz will be the requi.LU antenna isolation. 

The effective noise bandwidth B^ of the TACAN signal can be found 
from numerical integration of the measured spectrum, as shown in Figure 28. 
Thus, 


B 


0 



^ y S. Af = 0.7 MHz, 
^0 i=l ’ 


and the effective spectral density of the 1 kW (+60 dBm) TACAN pulse is 
60 dBm - 10 log 0.7x10^ = +1.5 dBm/Hz. Since the dynamic range of most 
spectrum analyzers is on the order of 60 dB and since it is necessary to 
determine power density levels approximately 180 dB below the TACAN den- 
sity level, extrapolation techniques must be used to determine the TACAN 
spectral density at fp + 77 MHz iGPS L2 center frequency). 

The first step is to ''ompare the measured TACAN power spectrum with 
the TACAN specification. According to the specification, the power in a 
0.5 MHz bandwidth centered tO.8 MHz from fQ (center frequency) and ±2 MHz 
from fQ shall be 23 dB and 38 dB, respectively, below the power in a 
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0.5 MHz bandwidth centered at fQ. By using numerical integration of the 
measured TACAN power spectrum, the power at ±0.8 MHz is found to be 46 dB 
below that at fQ, and the power at ±2 MHz is at least 60 dB (analyzer 
limit) below that at fg. Thus the actual performance of the TACAN, in 
terms of spectral shape, is at least 22 dB better than the specification. 
The problem now is to extrapolate this performance to Tq + 77 MHz. This 
is accomplished by utilizing time domain information. Examination of the 
meast»*ed time response of the TACAN pulse, as shown in Figure 28, indi- 
cates that the actual pulse shape is between a Gaussian pulse and a 
raised cosine pulse. In terms of spectral response, the raised cosine 
pulse is the worst-case interference source. The spectrum of the raised 
cosine pulse is given by 


r\a)/ -* VT p— 

WT [1 + (ut) ] 

3 

so that the power envelope falls off as the square of oj , or 

( 2 ). 

Thus, since the power in 0.5 MHz at fg is found by numerical integration 
to be 58.7 dBm, the power density in 0.5 MHz at 77 MHz, by extrapolation, 
is found to be 

P = 58.7 dBm - 60 - 60 log (-2-) - 10 log (0.5x10®) 

= -153 dBm/Hz. 

The required antenna isolation is thus 

Ant Iso = 181 - 153 = 28 dB. 

Several points should be noted. First, this analysis does not treat 
spurs which the TACAN transmitter might generate, since Axioniatix has 
not received a diagram of the TACAN frequency generator to analyze the 
potential spur emission problem. Second, this analysis has concentrated 
on the L2 GPS frequency. Since LI is 450 MHz away from the upper TACAN 
frequency, the wideband modulation sideband interference will not affect 
the LI GPS performance from a signal -to-noise standpoint. However, the 
GPS preamplifier must be .esigned so that the high peak power TACAN pulses 
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do not create serious intermoduiations due to saturation of the 
preamplifier. Thus, if the 28 dB antenna isolation presents a serious 
implementation problem, Axiomatix recommends ignoring the requirement 
and not using the L2 frequency for ionospheric correction when the TACAN 
is operating in the upper portion of the band. Furthermore, ongoing 
studies may find that the ionospheric delay correction is unnecessary in 
the context of the Shuttle navigation accuracy requirements. If this is 
the case, the TACAN wideband interference to L2 is a moot point. 

A description of the TACAN measurements made by JSC is given in 
Appendix A. 
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6.0 GPS PN CODE ACQUISITION 

The GPS signal structure utilizes an acquisition code, the C/A code, 
to enable signal acquisition to occur in a reasonably short time. None- 
theless, since the acquisition code itself is a PN code, of length 1023 
chips, the Shuttle GPS receiver must utilize a code search procedure to 
find the correct phase of the code for code acquisition. Two candidate 
code search and acquisition procedures for the Shuttle GPS receiver are 
discussed and analyzed here. The second of these procedures, the sequen- 
tial search technique, is found to offer a performance advantage over the 
first technique, the fixed step search. Consequently, the performance of 
the sequential search technique is optimized by optimum selection of the 
detector parameters. Furthermore, the effect of the Shuttle doppler on 
the search time is analyzed. 
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6.1 PN Acquisition Time Performance of the Shuttle/6PS Receiver in the 

Presence of Shuttle Signal Dynamics With a Fixed Step Size Search 

6.1.1 Introduction 

The Shuttle/6PS receiver is responsible for extracting the range 
and range rate data from the GPS transmitted signals. This is accomplished 
by cross-correlating the received signal plus noise with locally generated 
reference codes and searching for the position of maximum correlation. 

Upon receiving an indication of this maximum correlation position, the 
receiver then enters a search/lock strategy during which time a code 
tracking loop is enabled to complete the acquisition process and begin 
tracking the code. After removing the C/A and P codes from the received 
signal, the two L-band carriers can then be acquired and tracked by some 
form of suppressed carrier tracking loop (e.g., an AFC aided Costas loop). 
Finally, when both carrier and codes have been acquired and the corre- 
sponding loops are tracking, the data may be detected. 

The complete acquisition/track sequence (code and carrier) for each 
satellite on each L-band carrier must be performed in the presence of 
various Shuttle signal dynamics (e.g., center of gravity (C6) dynamics 
and dynamics about the pitch axis). The first question that must be 
answered is whether or not aiding is required and, if so, how much. 

This report begins an attempt to specify the required aiding by 
focusing on the code acquisition aspect of the acquisition/track problem. 

In particular, the PN code acquisition process is modeled as a single 
dwell time system as in [2], and expressions for the mean time to acquire 
code phase in the presence of doppler, doppler rate, and finite duration 
jerk (doppler acceleration) are developed. A single dwell time system is 
one in which only a single fixed time integration (dwell time) is utilized 
in the code search process, as opposed to a double dwell time search [2] 
which has two integration periods— one to quickly search for the correct 
code phase and the second to provide an improved estimate of whether the 
correct code phase has been found. The next section briefly reviews the 
Markov chain model of the single dwell time system developed in great 
detail in [2]. 
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6.1.2 Markov Chain Acquisition Model of the Single Dwell Time System 

Consider the simple model of a single dwell time PN acquisition 
system illustrated in Figure 29. Since the PN code is acquired without 
knowledge of the carrier phase, the model employs a standard type, non- 
coherent (square-law) detector. Briefly, the received signal-plus-noise 
is amplified and then correlated with a local replica of the PN code, 
the correlation being performed by a multiplier followed by a bandpass 
filter. The filter output is then square-law detected with the detector 
output being integrated for a fixed time duration (the "dwell time") 
in an integrate-and-dumo circuit. 

The Markov chain nature of the acquisition model stems from the way 
in which the integrate-and-dump (I&D) output is processed. In particular, 
if the I&D output is above a preset threshold, then a "hit" is declared 
and the system enters verification mode . If this hit represents a true 
hit (i.e., the correct code phase has been determined), then the system 
has officially acquired and the search comes to an end. If the hit is a 
false alarm, then verification cannot be consummated and the search must | 

continue. In either case, we shall assume (as in [2]) that the verification 
is characterized by an extended dwell time (e.g., kx^ sec; k>>l) and an ; 

entering into the code tracking loop mode. Understanding that a true hit 
corresponds to a single code phase position and this can only occur once 
per search through the code, we can regard the time interval kx^ sec as 
the "penalty" of obtaining a false alarm, since a false alarm can occur ■' 

on any code phase position. If the I&D output falls below the preset 
threshold, then the local PN code generator steps to its next position |- 

and the search proceeds. Thus, at each test position (aside from the 
single true code phase position), one of two events can occur, each char- ‘ 

acterized by a probability of occurrence; namely, a false alarm can occur 
with probability Pp^, causing a penalty of kx^ sec, or no false alarm 
occurs with probability (1-Pp^), resulting in only a single dwell time 
of Xp sec— hence, the Markov chain model. | 

The time to acquire T, i.e., the time to declare a true hit, is a I 

i ^ 

random variable and, in general, depends on the initial (at the beginning 

of the search) code phase position of the local PN generator relative to 

that of the received code. Since the probability distribution of T is !; 

difficult to obtain in practice, one ordinarily is content with determining 
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the mean acquisition time f and. In the absence of a priori knowledge 
concerning the relative code phase positions of the received and locally 
generated codes, the local PN generator Is assumed to start the search 
at any code phase position with equal probability. Stated In mathematical 
terms, letting q denote the number of cells (phase update positions of 
the local code replica) to be searched, then the probability of having 
the signal present (true hit) In the first cell searched Is 1/q, and the 
probability of It not being present there Is 1 - 1/q. For example. If the 
number of code chips to be searched Is denoted by N and the search proceeds 
In half-chip Increments, then q*2N and P^ = 1/2N. More generally. If It 
has been determined that the signal Is not present In the first k-1 cells, 
then the a priori probability Pj^ of finding It present In the j^th cell Is 
l/(q+ 1 - k), where (q+ 1 - k) Is obviously the number of remaining cells 
to be searched, each possessing an equal probability of having the signal 
present. 


6.1.3 Mean Acquisition Time Performance 


a. No Code Doppler 


Using the above Markov chain model, tne mean acquisition time of 
the single dwell system Is determined In [2], with the result (for 
large q): 


(2-Pjj)(l.KPp^) 


tp: 


q T 


D’ 


(i: 


where P^ Is the probability of detection (determining a true hit whea. 

In fact, the local code replica is at the correct code phase position). 
Furthermore, the result In (1) assumes that no code doppler Is present 
and that P^ Is time Invariant. Although Idealized insofar as Its assump 
tions, the result in (1) can be modified to reflect actual hardware or 
software systems. These modifications are discussed in [2] and are also 
germane to the extensions which follow. 

b. Code Doppler 

When code doppler Is present, the acquisition time performance 
of the system of Figure 29 is affected In two ways. First, the code 
doppler causes the relative code phase between received and locally 
generated PN codes to be time varying during the dwell time of the 
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Integrate-and-dump. This has the effect of Increasing or decreasing the 
probability of detection P^. depending on whether the code doppler Is 
causing the relative code phase to Increase or decrease. The second and 
potentially more dominant effect Is that code doppler affects the average 
search rate. In fact. If the code phase shift caused by the doppler over 
a single dwell time Is equal to the step size (phase update) of the 
search, then the average search rate Is reduced to zero. 

To take both of these effects Into account when computing mean 
acquisition time Is Indeed a difficult. If not Impossible, analytical 
task. However, It Is possible to account for just the effect of doppler 
on average search rate In a way which represents an obvious extension of 
the result In (1). Letting Af denote the code doppler In chips/sec and 
T^ the chip Interval In sec, then the mean search (code phase) update y 
In chips Is given by [2]; 

p - '^^c ^^c ^FA * 


where aT /T represents the search update In the absence of doppler or, 
equivalently, the step size of the search In fractions of a chip# 

Is the code phase shift due to doppler during the dwell time; and Afj.KTjj 
Is the code phase shift during verification caused by a false alarm. 

Thus, replacing q by N/y In (1) gives an expression for mean acquisition 
time In the presence of code doppler, namely L2], 


T 


2Pr 


(2..Pp)(UKPp^)NTQ 

'^♦if^Xo(l*KP„r 


(3) 


Since can be either positive or negative, then depending upon Its 
sign, the code doppler can either speed up or slow down the search. 

With regard to the magnitude of the code doppler, we shall assume that 
|Pfcho(l*KPp ) << aTc/Tc, so that the denominator of (3) never approaches 
zero and the search always proceeds in the direction dictated by the code 
phase update provided by the local PN generator. Finally, note that, when 
Af * 0, then aT /T = N/q and (3) reduces to (1) as It should. 
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c. Code Doppler and Doppler Rate 


The way In which code doppler is accounted for in extending the 
result of (1) to that given in (3) can also be applied to further extend 

(1) to include the effect of code doppler rate. In particular, we shall 
compute the mean search update analogous to (2) and again replace q by 
N/u in the expression for mean acquisition time. 

When code doppler alone was present, we observed that the mean 
(statistical) search update was time invariant, i.e., the expression in 

(2) characterizes every cell being searched. When, in addition, code 
doppler is present, the mean search update is now time dependent in the 
sense that it is now a function of the cell being searched. Thus, we 
start our :malysis of this case by examining the mean search update in 
the first and second cells from which we shall find a general expression 
that characterizes this quantity in the ith cell, 1 <.i <.q. 

In general, the mean search update u in any given search cell is 
equal to the nominal search update AT /T (typically, 1/2 for half-chip 
search increments) of the local PN code generator plus the mean change in 
phase of the received code over the search time of that cell. Thus, 

letting p, denote the mean search update iii the first cell being searched, 

. ‘ 2 
ai.d Af the code doppler rate in chips/sec , then from the Markov model 

previously established for the single dwell time system, we have that 


if * - "fa' 


rtf 


dt + Afc Pp^ 


r(K+l)' 


dt 


^ ^ ^ ^^c ^FA I 


(K+1 )t. 


t dt 


AT 


T ^ ^^c ^ ^FA^ * 2 '^^c ^FA^ + 2K) ] . 


(4) 


Similarly, in the second cell tr be searched, the mean search update 
would be given by 




♦ 

We continue to assume, as before, that the code phase derivatives 
are positive when they are in such a direction as to aid the search 
(reduce the acquisition time). 
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AT r f2r^ f(K+2)TQ 

T7*‘V'-'’fa> ('-'■fa) , ■“"'■faL 

w L Tn ip _J 


* ''fc"FA 


♦‘ V'-W 


f(K+2)x 

f2(K+l)T 

(1-Pp.) dt + Pp. 

MK+l)Xp 

° dt 

MK+l)Tp J 


(1 


r2x f(Kt2)'n 

-^ a )I/““*' faJ , ““ 

D ~ - 


(K+2)t, 
0 




• ^^"Af,TpO*KPFA) + ^if<,Tp2[l*P„<''^*2K)] 


(’-^ a'I 


0 ,2(Ktl)tp 

t dt + Pp^ ' ^ 


(K+2)t 

t dt + Pp. I 

(K+1)t„ *^^J(K+1) 


f2(K+l)T 

^ t dt 

J(K+1)t. 


+ A^J. [Tp (1 + KPp^)] . 


(5) 


To generalize (4) and (5) to an arbitrary search cell, we first observe 
that, for any time tp > 0, 


rt«+rp |.tp+(K+l)Tp 

(1-Pfa) 1“ °1t*Pp^ [ “ = Tp(fKP„) 


( 6 ) 


and 






t dt + P 


FA 


to+(K+l)Tp 


tdt 




= tTp[l+KPp^] + iTp2[l + Pp^(K^ + 2K)] 


(7) 


Thus, extending; the Mar'^ov chain procedure leading to the expressi r. for 
^2 and using (6) and (/), we arrive at the general expression for the 
mean search update of the nj^st cell, namely, 

VI " 'D <' * '"'fa' * 7“^ 'D t' * ''fa 

C 

♦ Af, ,p (U KP„) I Pf\(l - Pp^)"-'‘ [(n - k) * k (K* 1 )] (8) 
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y 
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t 


? 


1 


* 


I 




> 


where 


r ^ nl 
rn ■ k! (n-k)! 


(9) 


1s the binomial coefficient. The third term In (8) can be further ampli- 
fied as follows: 


L k<:n "fV - '■fa)"'" C" • J„k'n'’FV’-'’FA'""'("*'‘''i 


n-k 


k»0 


k=0 


Since, from the properties of a binomial distribution, 

k<:/FA - '’fa'"'' ■ ’ 


( 10 ) 


^ ‘ np,,. 


k=0 


k'n^FA^' '^FA 


FA 


( 11 ) 


then from (10) and (11), we have 

I * n + KnPp^ = r.(l + KPp^) . (12) 

k— 0 


Finally, substituting (12) Into (8) gives the desired simplified result: 
“ntl “ >D<' * K ’’fa' * ^0 [’ * ’’faI"^ * 2K)] 




(13) 


Note from (13) that the mean search update is a linear function 
of the search cell. Because of this dependence on n , we cannot directly 
replace q by N/p in (1) to arrive at a formula for mean acquisition in 
the presence of code doppler and doppler rate. Rather, we shall first 
find the average mean search update u obtained by averaging p of (13) 
overall q search cells, i.e.. 



1 l-J 
*1 nio 


(14) 


and then make the above suggested replacement in (1). Thus, from (13) 
and (14), we have that 
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= ^ * ‘fc 'D<' * '< ’’fA> * i 'D t' * ’’fa 
*( i^)Af,[xD(l*KP„)]' 

' '''cV'^'’fA> I<'('*'<'’fA>^*''^'‘fa(’-'’fA>I- <’=) 

Since the total mean search update must correspond to the total number of 
chips searched, we have that* 

‘1:"' 

= I u , 1 = N. (16) 

n=0 " ‘ 

Thus, substituting N/a for q in (15) results in a quadratic equation 
in M, namely, 

■ ’ Ilf' ^ *’fA^ * I ^^c ■ ^FA^ 

- ]-AfcTp^N(l+KPp^)^ " 0. (17) 

Letting 

A = 1 

® f “ ^^c W ■*■ 2 ^^c ^FA^^ " ^FA^ 

C = 1 Af^^P^N(l+KPp^)^, (18) 

then 



here we have made the assumption chat jAf^l anu |z\fj.j are small 
enough such that un+i > 0 for all n = 0,1 , . . . ,q-1 and thus the search pro- 
ceeds only in one direction, namely, that dictated by the local PN gener- 
ator code phase update. 


I 


Finally, substituting N/u for q in (1), with y defined by (19), gives 
the resulting expression for mean acquisition t%e in the presence of 
code doppler and code doppler rate, namely. 


T = 


(2-Pd)0*KP„)N,o 


2P, 


[uJTTf) 


2A 


( 20 ) 


Note again that, when Af = Af = 0, we have that B = aT /T = N/q, C=0, 

c c c c 

and (20) reduces to (1). 


d. Code Doppler, Doppler Rate, and Doppler Acceleration 
O ver a Fini te Time Interval 

When, in addition to code doppler and doppler rate, the received 
signal dynamics are characterized by a doppler acceleration or jerk over 
a finite interval of time, say, 0 > t •: t^ , then the received code phase 
e/2;i takes the fonn 


0 1**21"' 
^ + Af t + -^ Af t^ + ^ Af t‘ 
2ti c 2 c 6 c 


0 . t V t^ 


2T 


+ Af t + 4 Af t^" + i Af t' 
2 it C 2 c 6 C 


'6 ■ ^1 ^ ’ ^1""^ 


( 21 ) 


where Af^, denotes the doppler acceleration. The second equation in (21) 
ca- written in a more convenient form which identifies the effective 
code phase, doppler, and doppler rates, namely. 


S ' S * 5 ''fc ‘ ' * <-''c - 2 ‘ • 't - ' 


J- + Af ' t + 1 ' t\ 

2n L 2 c 


( 22 ) 


For convenience jf analysis, we shall assume that the time dura- 
tion of the jerk is an integer multiple of the df'"'! time, i.e., t^ 

M 1, and furthennore, less thari the time duration penalty of a false 


alarm, i.e., < (KH)xp or, equivalently, M<K+1. As before, we shall 

characterize, for the purpose of illustration, the mean search update 
in the first and second search intervals from which an expression valid 
for the Uh search interval will become obvious. Finally, the (time) 
average mean search update will be computed and used to evalute the mean 
acquisition time performance of the single dwell system in the presence 
of code doppler, doppler rate, and finite time duration jerk. 

In the first code search interval, the mean search update is 
now given by 


Tj. ^ 


rT 


dt + Af Pp. f dt 
c FA Jj 


(K+1) 

D 

”"n 

^ ^t dt . Af^' Pp, J 


^ ^^c PpA 


" ^^c '^FA 


f(K+l)TQ fTrj 

°dt + Af (1 -P .) ^t dt 
Mt,, ^ Jo 

r(K+l).o 


tdt 


'Mt, 


Mx 


FA 


^t2 dt 


AF 


+ Af, Tn(l+KPp^) + ^ Af^ip^n +Pp^(K^ + 2K)] 


■ “'c ‘D 


■ 2 ^^c” ~ ^FA^ ^FA^^} 

+ 2 Af^ M Tp |[1 + Pp^(K + 2K)] - [(1 - Pp^) + Pp^M 


* 6 "^FA^ ^Fa”^^ 


(23) 


Note that, when Af^ = 0, then (23) reduces to (4). In the second search 
interval, the mean search update is similarly given by 



f 






t 


♦ 


» 
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-rf^n-Pp,) 




2x 


"0 


'^(Af^+Af^t+ i Af^t^) dt 


+ P 


FA 


Mt 




(K+l)t 


^(Af +Af t+ I Af, t^) dt 
' c c 2c' 


Mt 


° (Af ' + Af 't) dt 
C C 


+ P 


FA 




D 

(K+2)t 

^ (Af ' + Af 't) dt 


+ p 


FA 


2(K+l)i 

(K+l)xr 


(Af ; + Af 't) dt 


(24) 


which upon simplification becomes 

* Af,. t D<' * ‘2 ''^c*''D 

« {(l.KPf^)-(l-Pp,)[(l-P„)*Pp,(M-1)]} 

+ tKPFA>^t j['*P„(K^+2K)]} 

-IaFcM.^o-Pfa) [30-Pp,)*Ppy-in 

"1''c'0('-'’fA>[^('-'’fa'*'’fa‘'’^-’»- (25) 

Once again, when Af^=0, (25) reduces to (5). Extending this Markov 
chain approach to the n+l st cell, we arrive at the general expression 
characterizing the mean search update n^^-j of this cell, namely. 


I 
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c 


n/^fcC^oO + KPp^)]' - + KPp,) - (1 - Pp^)"[(l - Pp^) ^ PpA(M-n)]} ; 


+ Af^Mx3|n{l + KPpA)2+ }[1 +Pf^(k2+2K)]} 

- ■ P«)"[0 -PpA)(2r> +1) +Ppa(M^ - n^)] 


+ -Ppj"[(l-PpJ(3n^ + 3n+l) + Ppfl(M^-n^)]; 0<n<M-l 


FA' 


FA 


FA' 


(l-Pp;^)(n+l)^+Pp^M^-n^ 


(26a) 


and 


'n+1 


' ^ " ‘V ^dC * "'’fa' " J *V’D [' * I’fa'"'" 2K)] 

c 


+ HA^^' [Tp(l +KPp^)]^ 


n > M . 


(26b) 


Using the definition of average mean search update y given in 
(14), then substituting (26) into this equation gives 


aT 


p -■ 


^ * A V T„0 * KPp^) * 1 ,2 D * Pp,(k2 . 2K)] 


“a q-1 “ 
q ^ " 

^ n=0 




» [0-P„)*HPfa] 


~i M-1 

? ^ <’-^A> 
^ n=0 


- P 


FA 


1 M-1 r 

I I nO-%)" 


n=0 


lAfcMx3|[l.Pp,).M-Pp,] 

Pi n 

■*■ ^(^ ■ PpA^ n I " PpA^ 


n^O 
M 


1 

^ n=0 
P 


I (1-Pp/ 


FA 


1 M-1 9 


311 V 3/, D 1 V n3/i D ^'14 .m3d 
6 J q L ” ■ ^FA^ ‘ q _L” “ ^FA^ ^ ^FA 


n-1 


FA' q 


n=0 


1 nil 


(27) 



Froin the summation properties of arithmetic and geometric series, we 
have that 


q-1 

n=0 


2 


M-1 

I (1-Pfa) 

n=0 


’’fa 


Yn(l-P )" = 


FA^ 


FA' 


n=0 

M-1 


FA 


FA 




l n (1-W" “ 


n=0 


FA' 


- (M-l)(l -Pp^) 




FA 


n(i-Pf/-'* (i-Pp/l 


FA 


I r^d-Ppj)" - V »^0 -Pfa)" = m3(1-P„)". 

n=l n=0 


FA' 


(28) 


Substituting (28) into (27) and simplifying yields 
'c 


li = ^0^^ ■^'^*’fA^‘^ 2^V^D ^FA^^ ■ ^Fa4 


, 1 * £ iui3 3 

* ’o ■ 


(29) 


Comparing (29) with (15), we observe that the effect of the finite dura- 
tion doppler jerk on the average mean search update is simply the addi- 
tion of the term Af^M Tp/(6q), which is inversely proportional to the 
number of search intervals q . 

Letting q=N/p as in (16), we once again get a quadratic equation 
in y, namely. 
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I 


» 


I 


> 


t ^ ^ * ''^C 'o<' * ^"fa' ^ 1 ‘V -0 ’’faC - '■fa> 


- 0 . 


(30) 


Letting 


A = 


B = 


C = 




AT, 


jf * At,’ T„ (1 * K Pp,) + I A^, 2 P„(l - Pp^) 
N()*KPp/. 


(31) 


wi 


1th, as a reminder [see (22) with t^ = MtQ] 

^ Af m2 2 

J - 2 Afc" ^0 


if = if. - i if. ^ 


ifj,' = if^tif^HTp, 


(32) 


then the mean acquisition time f is once again given by (20) with A, 

B, and C now as defined in (31). 

6.1.4 Evaluation of Detection Probability Pq and False Alarm Proba- 
bility Pp/\ in Terms of PN Acquisition System Parameters 

The formulas for mean acquisition ti.ne developed in the previous 
section are all functions of the detection probability Pp, false alarm 
probability Pp^, and dwell time t^. Thus, it would appear at first 
glance that, for specified values of detection probability and false 
alarm probability, one could arbitrarily select the dwell time to 
achieve any desired mean acquisition time. Upon closer examination, 
one realizes that ndeed this is not possible with the fallacy lying 
in the fact that, for a given Pp^ and predetection signal -to-noise 
ration, P^ is implicitly a function of t^. To place this statement 
in evidence, we begin by evaluating Pp and Pp^ in terms of the PN 
acquisition system parameters for the simple case of no code doppler 
or doppler derivatives. Following this, we shall generalize these 
results to correspond more realistically to the Shuttle signal dynamics. 
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a. No Code Doppler 

When signal is present (i.e., the cell being searched corresponds 
to a sample value on the PN correlation curve), then the input to the 
square-law detector can be expressed in the fcrm 


x(t) = s(t) + n(t) 


where 


^ A cos (o)Qt+ til) + *^ n^(t) cos (o)Qt+ i|i) 
- ^ n^(t) sin (u)Qt + tii) 

^ R(t) cos (wQt + t|i+ e(t)) , 


R(t) = J(A + n^(t))^ + ng^(t) ; e(t) = tan 


-1 


ng(t) 

A+ n^(t) 


(33) 


(34) 


★ 

In (33), A is the rms signal amplitude, wq the radian carrier frequency, 

and n (t), n^(t) are bandlimited, independent, low pass, zero mean 
c s p 

Gaussian noise processes with variance o^^NqB, where Nq is the single- 
sided noise spectral density and B Is the noise bandwidth of the pre- 
detection bandpass filter. 

The output of the square-law device in Figure 1 in response to 
the input x(t) of (33) is (ignoring second harmonics of the carrier): 


y(t) 4 x^(t) = R^(t) = (A+n^(t))2+n2(t) (35) 

and has a first order probability density function (pdf) which is given 
by 


p(y) 



0 ; otherwise 


y >. 0 


(36) 


For the moment, we shall not enter into a discussion concerning 
the various system losses and gains which ente^* into the calculation of 
the effective signal amplitude to be used in predicting true system 
signal-to-noise ratio behavior. Such a discussion will be given later 
on in the development. 



where IqU) the zero order modified Bessel function of the first 
kind and 

Y A aL 

^ NqB 


( 37 ) 


is the predetection signal-to-noise ratio. In the absence of signal, 
i.e., A=0, (36) reduces to the chi-squared pdf: 


p(y) = 


1 

2 

0 

0 ; 


exp 




otherwise 


y C 


(38) 


which characterizes the square-law output in all search cells which con- 
tain noise only. 

If y(t) is sampled at intervals T= 1/B, then these samples are 
approximately independent, and furthermore, the Integra te-and-dump 
output can be approximated by a summation over these sampled values, 
namely. 


2 4 JL 

I 


rTr 


Ng-l 


y(t)dt = ^ I y(kT) 
”b k=0 


(39) 


where* 


"b - T = 


(40) 


Using the approximation in (39) and the first order pdf's of 
(36) and (38), the pdf of Z, namely, p(Z), for signal present is 
given by 


P(Z) = 


\ /J- 

a2 1,0^ 


v(Nn-l)/2 

— — 
/ z \ 


2 Z 

) exp 

_■ b(„2 ■ V 

'Nj-I 



Z > 0 


(41) 


0 ; otherwise 
and for signal absent is given by 


★ 

It is convenient here to assume that Ng is integer, although the 
results which follow are, for large Ng, valid for Ng noninteger. 
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P(Z) = 



< 


Z 0 


(42) 


0 ; otherwise 

In (41), Ing.i(x) is the (Ng-l)st order modified Bessel function of 
the first kind. 

2 

Normalizing Z by o /Ng=NQ/Tjj or, equivalently, letting 
Z*-ZNg/a^, we can rewrite (41) and (42), respectively, in the 
simpler form. 


P(Z*) = 


/ 1* r -1 

(ib) «p (- z* - "b v) In,-i [zfB ^ z*J : z* ^ 0 


V 

0 ; otherwise 


(43) 


and 


p(z*) = 


( Nb -')1 
0 ; otherwise 


exp (- Z*) ; Z* > 0 


(44) 


The probability of false alarm, Pp^, is the probability that Z 
exceeds the threshold n when signal is absent or, equivalently, in 
terms of the normalized random variable Z* and the normalized threshold 
n* - nNg/o^, 


FA 


p(Z*) dZ* = 1 


= e 


. n * , ,( Nb -‘')/2 

0 (Ng-1)! 

Nr- 1 k 
1* y (n*) 

k=0 


exp (- Z*) dZ* 


k! 


(45) 


The detection probability Pg is the probability that Z exceeds the 
threshold n when signal is present. Thus, using (43) rather than 
(44), we get 
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it i V I 

0 \V/ 


(Ng-l)/2 


exp(-Z*-NgY) , [2^NgYZ*JdZ* 


( 46 ) 


which cannot be put into a closed-fonn expression. 

For large Ng (the case of most practical interest), things become 
quite a bit simpler. Defining Yl* * y(kT)/a^ , then from (36) and (38), 
the pdf's of yj^ in the presence and absence of signal are, respectively, 

r -(y|^+Y) 


p(yL*) = { 




y.*> 0 


?(y,*) = 


0 ; otherwise 


; y^ ^ 0 


(47) 


0 ; 


otherwise 


(48) 


Also, from (39) and the definition of Z* in terms of Z, we have 

Ng-1 

z* = I y, ‘ (49) 

k=0 

Since, by previous assumption, the y|^*'s are independent random vari- 
ables, then for large Nn, Z* is approximately Gaussian distributed with 

° 2 2 

mean Z* = Ngy* and variance = Ngo^^ . The means and variance of 

the pdf's in (47) and (48) are well known [3] to be 


/* = 


/* = 


1 + Y ■» 0 * = 1 + 2y ; signal present 

y 


= 1 


signal absent 


(50) 

(51) 


Thus, 


I* = 

I* = N 


Ng(l + y); <^ 2 * ^ Ng(l + 2Y) ; signal present (52) 

2 


B 


°Z* " 


; signal absent . 


(53) 
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Using the Gaussian assumption, the false alarm probability is 


FA 


r 1 

■ (Z*-Nb)^’ 

^ n* Ng 

2Ng 


dZ* 


= ^ erfc 


n* - Ng 


/71 


B-I 


where 


j erfc $ , 


(54) 


erfc X - 


_ 2 _ 


[ exp (-t^) dt. 


(55) 


Thus, if Pp^ is specified, 6 can be determined. The corresponding 
detection probability under the same assumption is 



II 

1 

a. 

r - - exD 


f 

J 

*n* y^irNgd +2 y) 

2Ng (1 + 2 y) 


dZ* 


t 


» 


= j erfc 


n*-Ng(l + Y) 


2 erfc 



(56) 


» 


I 


Combining (54) and (56) and reidentifying Ng and y in terms o 
system parameters gives the final relation 



Thus, given Pg, 



a id B, the dwell time 


^0 


is determined. 


t 



Before using (57) and the dwell time determined from it in the 
formulas derived in the previous section for mean acquisition time, 
several modifications based upon practical considerations must be 
made. 

(1) E ffective Probability of Detection and Timing Misalignment 

The calculation of detection probability as in (46) or (56) 
implicitly assumed that only one cell in the entire search satisfies 
^he "signal present" hypothesis. In actuality, since the PN correla- 
tion curve exists over an interval of ±1 chip around the peak, then, 
for example, a system whicn updau the locally generated code phase 
in half-chip increments would yield several cells for which signal 
could be considered present. 

Typically, the system is designed on the basis of the worst case 
correlation, which for the half-chip jpdate case would correspond to 
the pair of correlation points ohe-quarter chip away from the correla- 
tion peak. Since the normalized correlation value at these points is 
0.75 (relative to a peak of 1), then the single signal point calcula- 
tion of detection probability as in (46) or (56) would be based on an 

2 

effective reduction in the nominal signal-to-noise ratio A /Nq of 
10 log^Q(.75)^ - 2.5 dB. Since, however, in reality there exist two 
worst case correlation positions, then the effective probabili ty of 
detection Pp for use in computing mean acquisition time is computed as: 

’’o' = "d* = 2'’d- 'd • 

where the first term in (58) represents the probability of detecting 
signal present on the first, correlation pcint and thp second term is 
the joint probability of not detecting signal on the first correlation 
point and detecting signal p esent on the second correlation point. 
Clearly, for low signal-to-noise ratios (small P^), the effective 
detection probability is approximately twice that computed on the 
basis of a single signal present cell [(46) or (56)]. In summary, 
then, the computation procedures would be as follows. For a given 
Pp^, determine 3 from (54). For a specified P^ , find Pp from (57), 
degrade the given nominal value of y by 2.5 dB and solve for Ng in 
(56). Determine the dwell time f“om Tg=Mg/B, wht'". B is the given 



bandpass filter bandwidth dt,_rm1iied by considerations on allowable 
modulation distortion (to be discussed next). Using , Pp^, and Xp 
In (1), solve for T. 

(2) Mod 'J lation Distortion Effects 

Typically, the PN modulated carr^T Is also biphase modulated 
by data. Depending on the ratio of predetection filter bandwidth B 
to data rate R, this data modul.it^of’ will suffer distortion and an 
equivalent power reduction as It passes through this filter. The 
equivalent power reduction factor D is computed from 


D 


S^; *) |H(j2irf)|^df . 

Jm 00 


(59) 


where Sjjj(f) Is thv power spectral density of the data "Podulatlon and 

H(j 2 tt f) Is the equivalent lowpass transfer function of the precietec- 

2 

tion bandpass filter. Thus, the nominal signal power A must be multi 
plied by D to account for this effect when computing the effective 
s1gn'l-to-noise rat-o to be used in xhe previous detection and false 
alarm probability computations. 


(3) Reduction in Noise Spectral Density Caused by PN Despreading 

Multiplicat'on of the equivalent nolee process at the PN acqui- 
sition system Input by the locally generated PN sequence spreaos the 
spectrum of this noise process and simultaneously reduces Its effective 
spectral height Into the data filter. Letting Nq denote this effec- 
tive noise spectral density, then since the bandwidth of the data 
f predetection) filter is much narrower than that of the PN process, 
we have that 


N 


0 


(•“ /sin 7t f T 



(60) 


In (60), we have assumed for simplicity of the calculation that the PN 
line scectrum is approximated by its envelope. Thus, again in computing 
the effective signal-to-noise ratio of the system, should be replaced 
by Nq' . 

Finally, summing up the effects of 6.1 .4.a(l ) - 6.1 .4.a(3) , the 
effective signal-to-noise ratio y' in the predetection filter bandwidth 
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t 


is given by 


y‘ 


A^DL 

upr* 


(61) 


where D is defined in (59), Nq' in (60) and L, the loss due to a frac- 
tional chip misalignment A'T^ from the correlation curve peak, is given 
by 



Age n, for I^euf-chip search updates, the worst case loss corresponds 
to .VTe-'/4. 

b. Code Doppler a id Its Derivative s 

When code doppler and its derivatives characterize the received 
signal dynamics, then the timing error between received and local PN 
codes is not constant over the dwell time of the detection process. 

This "smearing" effect has a direct bearing on the calculation of the 
loss due to fractional chip timing misalignment and thus the expression 
for this loss given in (62) requires modification. 

Consider first the case where the received signal is character- 
ized by code doppler and doppler rate. Then the normalized signal input 
to the acquisition dwell time integrator in Figure 29 is (under the signal 
present hypothesis): 

s(t) = ^1 - -y- - Af^t - ^Af^ t^^ . (63) 


The corresponding normalized dwell time integrator output is 


L = ^ 


0 


s(t) dt 


A’ T 


1 


A' T. 


“ Tq j ~ 'D 3 'D ^ y ■ T^ 


+ 'd ^ ^ ^c ^ ^c ^ ^ ^c 
3 4 20 


( 64 ) 
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Por small Af and Af , the loss L of (57) is well approximated by 

\0 


L * 


1 - 




i2 


( 65 ) 


Finally, if the "signal present cell" is being searched during 
the finite time interval over which the doppler jerk is present, then 
(63) becomes 

\2 


A' T 


c A -p 4- *.2 1 (66) 


s(t) = ^1- -Af^t-^Af^t - g Af^t 


and the corresponding loss as computed from (64) is 


3 3 




where, for simplicity of notation, we have defined 

AfJ’^ = Af. 






(68) 


Again for small Af^, Af^, and Af^, (67) is well approximated by 


L -- 


A'TX 


r,2 


(69) 


Using (67) or (69) in (61), one is able to compute, to a first 
order approximation, the effective signal-to-noise ratio in the predetec- 
tion filter bandwidth when the received signal dynamics are character- 
ized by code doppler and its first two derivatives. 

Because of the dependence of L on dwell time t^, the procedure 
for ultimately calculating mean acquisition time is more complex but 



may be summarized as follows: 

(a) For a given detection probability (Pq) » solve for 
from (58). 

(b) For a specified bandwidth B and filter type, compute D 
from (58) and Nq from (60). (Actually, carrier and code doppler 
affect these computations but we shall assume these are second-order 
effects. ) 

(c) Using 0 and NJ computed in (b) and L from (67) or L from 
(69), determine y' f>^om (61) as a function of Tq (assume A^/N^ is 
given and worst case misalignment). 

(d) Letting Ng=BTp, then for a given Pp^, together with Pp 
determined from (a) and y' from (c), solve (transcendental ly) for 
from (56). 

(e) Using Tjj as determined in (d), and the given values of 
Pp and Pp^, calculate T from (20). 



6.2 Sequential Detection Acquistion Time Performance 

of the Shuttle GPS Receiver 

6.2.1 Introduction 

Sequential detection is a detection process first used in radar 
and, more recently, applied to the acquisition process in PN receivers. 

It offers the advantage of significantly faster search times when com- 
pared to "fixed step size" search procedures. The advantage is achieved 
because the sequential process is able to decide that signal (i.e., PN 
correlation or partial correlation) is absent much faster than the fixed 
integration time process. This is because the latter process spends as 
much time deciding there is no signal (which is the predominant test in 
an acquisition search) as it does in deciding that the signal is present. 
Thus, most GPS receivers are implemented with a sequential detection 
acquisition search algorithm. The following performance analysis of 
sequential detection acquisition is presented for the purpose of optimiz- 
ing the Shuttle GPS receiver acquisition algorithm parameters. This 
results in the minimum Shuttle GPS receiver acquisition time and, conse- 
quently, minimum time to first fix (TTFF). 

6.2.2 Description of Sequential Detector 

The sequential synchronization detector used for GPS receivers is 
shown in Figure 30. The mean output of this system is proportional to 
the total power contained in the input. The actual output X(t) is com- 
pared to a threshold e' that is chosen to maximize the probability of 
detection, P^, for a given probability of false alarm, Pp^. This will 
lead to minimization of the average synchronization time. 

The integration of the square-law output before making a threshold 
comparison increases the probability of detection, Pp, without changing 
the probability of false alarm, Pp^^. The bias subtraction before the 
integration process will cause the post-detection integration (PDI) out- 
put to increase when signal -plus-noise is present and decrease when noise 
alone is present. The threshold 0* is set to a negative value at the PDI 
output. If the output is above the threshold for all 0<.t<.Tj|^, where 
Tjp denotes the truncation time, the signal is declared present; other- 
wise, the output will go below the threshold, causing a dismissal of 
noise alone in a time much smaller than T^p, as shown in Figure 31. 
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PDI 

Output 


Signal + Noise 




Noise Only 


Figure 31. PDI Output as a Function of Time, t 


6.2.3 Signal Description 

The received GPS signal r(t) is of the form 

(‘^t + 4>) + d(t) (t) cos (o)t + <()) + n(t) 


where 

Pp(t) = high-precision ranging PN code of rate 10.23 Mc/s 

Pc/^(t) = Gold code of length 1023 chips and rate 1.023 Mc/s 

d(t) = 50 sps data signal 

0 ) = carrier frequency 

4> = unknown carrier phase 

P = received power with P^/Nrt=31.6 dB-Hz 
P P 0 

^C/A * •^eceived power with dB-Hz 

n{t) = additive thermal noise assumed to be Gaussian with 
flat spectral density kNg/2 . 

There is also a code doppler of ±5 chips/sec. 

6.2.4 Detector Description 

As shown in Figure 30, the receiver multiplies r(t) by a replica 
of p(t), which starts at some arbitrary initial code phase. The result- 
ing product is passed through a filter wide enough to pass the data 
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signal d(t) plus the doppler shift with little distortion. A filter 
bandwidth, B, of four times the data rate is enough. For a 50 bps data 
signal, a bandwidth of 200 Hz is chosen, although the effect of changing 
bandwidth on the results is discussed later in this report. If the 
initial phase of the code p(t) is very close to that of the received 
code, the filter output Z(t) will be the data modulated carrier plus 
noise; otherwise, the product will be noise alone. 

The noise, n(t), in (70) can be divided into its two orthogonal 
components, yielding: 

n(t) = nj(t) cos (ut + ij>) + ^ sin (wt + (<)) , 

where nj and ng are two independent, zero mean, Gaussian noise processes 
with spectral density Ng(in) = 1.18 Ng when in correlation and Ng(out) » 
1.4 Ng when not in correlation. 

Assume the square-law output takes the form v(t) = l/2 Z^(t). 
Then, the outputs under the two hypotheses are: 


under H 


0 ’ 


v(t) = I p2p d{t) + + ng^j-^^j(t) under . 

Let v(t) be sampled at t= 1/B; then the resulting random variable 
(v') will have the densities 


P(v'|in) = ^ I«(2 ^av' ) ; v'>0 under H 


1 


P(v' lout) 


= lp-v7R. 


V* > 0 


under H 


0 ’ 


where 


V* = 


Ng(in) B ' 


v' can be generated using the computer [4] as follows: 

.2 . .2 




, for the in condition ; 



v' s 2 » for the out condition ; 

" • 

and 92 are Independent, zero mean, unit variance random numbers and 
a Is the signal-to-nolse ratio In the filter bandwidth B. 

For the sequential detector to be effective, the bias (b) must 
take a value In between the means of the square-law output under the two 
hypotheses, I.e., 


E(X|out) < b < E(X|1n) . (71) 


where E[X] is the mean of the signal after bias subtraction and is given 
by [4]: 

E[X|out] = (1 - b’) NQ(out) B 
E[X|1n] = (1 +a) hQ(1n) B- b'NQ(out) B . 
Substituting in (71) yields: 


As a better value for simulation, take 


then. 


b = b' 


NQ(out) 


NQ(out) 

NQ(in) 


< b < 1 + a 


or 1.19<b<l+a. 

Although the optimum bias b can be theoretically bounded by the 
above values, the values of threshold 9 and truncation time cannot. 
However, adequate values for 6 and can be chosen based on simulation 
results such that the average synchronization time reaches a minimum 
value under a desired Pq, Pp^, and signal-to-nolse ratio a. Evaluation 
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of probability of detection Pp and average time to dismiss the noise alone 
Tq are needed to obtain complete information about the average synchroni- 
zation time T . This can be illustrated as follows, using the C/A code. 
Assume the length of the code 1023 chips and the search process proceeds 
in half-chip increments; then the number of cells to be searched (the 
uncertainty region) will be 2^ 1023 cells. The synchronizer takes an 
average of 2xl023xTp seconds to search through the whole uncertainty 
region. Since Pq is small for small SNR, it will take, on the average, 
1/Pp passes through this region before the correct synchronization point 
is reached. Thus, the average synchronization time is given by 


T 


SI 


2 X 1023 



(72) 


The values of T^ and P^ are calculated from computer simulations 
using the Monte Carlo method as follows. Let XN denote the total time 
taken to dismiss up to the truncation time T^j^ for a number of trials 
equal to XTEST, assuming the noise signal alone is present. Let Xj denote 
the number of trials that were not dismissed (the number of trials where 
the PDI ''utput is above threshold for all the time 0£til) under noise 
signal alone. Let Xjl denote the number of trials that were not dismissed 
under signal + noise conditions. Then, 

Tp = (XN)/(XTEST-Xj) 

Pp = (Xjl)/(XTEST) 

Pp;^ = (Xj)/(XTEST) . 


It is important to note that there is an average synchronization 
time correction (T^s2^^S3^ added to the average synchroni- 
zation time Tg^ given in (72). T^^ due to false alarm causing the 

test to proceed to truncation time and is given by 

T^2 = 2 X 1023 X T.j.j^ sec x sec . 


Tj 3 is due to time (equal to 
determine the false condition. Thus, T 


spent in the phase-lock loop to 
is given by 


T$3 = 2x 1023xTpg^g^^^x pp^ 
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where Is 1 sec, which Is the typical code loop response time, 

so the total average synchronization time would be Tg “ ^S1 ^^S2^^S3' 

total 

6.2.5 Simulation Results and Parameter Optimization 

A computer simulation of the sequential detector Is used here, 
since the sequential detection process Is not amenable to analysis. A 
flow diagram for the basic computer simulation program Is shown In 
Figure 32. The subloops for optimizing the values of threshold 6 and 
bias b are not shown, for purposes of clarity. For different values of 
threshold 6, the bias b Is changed within the permissible range so that 
a best value for b and 6 can be selected for fixed values of normalized 
truncation time ITRNK of 20 and number of trials of 1000. The relation 
between average synchronization time with b and e Is shown In Figure 33 
In the range where probability of detection (assuming single-signal 
level) Is not less than 0.6 and '■fa not greater than 0.005. The optimum 
values for b and e are chosen from curves to be b=2.4 and e = -1.2. Then 
a trial Is made to get the relation between truncation time ITRNK and 
average synchronization time T^-j. Two values of ITRNK are chosen to be 
16 and 10 where, as shown In Figure 34, ITRNK of 10 gives a smaller Tg^ 
but a bigger Pp^ than that of 16. Figures 35 and 36 show the relation 
between SNR(a) In dB versus Tg-| for the ITRNK of 10 and 16, respectively. 
As SNR In the receiver drops, the acquisition time Increases because the 
receiver nas to make several passes through the frOl period until synchro- 
nization Is detected. 

6.2.6 Correlation Peak Detection 

A trial using the actual signal which. In the Ideal case, takes 
the triangle shape shown In Figure 37 Instead of using single-point 
signal + loss was Investigated. As expected, the Pp Increases (e.g., 
for the previous specific values of b=2.4, 0»-1.2, ITRNK* 10, the 
signal point Pp= 0.686 while, for the correlation signal case, Pjj= 0.936). 
Table 19 shows In detail the actual Monte Carlo random starting positions 
and the associated probabilities. The Initial condition Is chosen at 
random and the trial Is repeated five times (L*5) with the resultant 
probability of detection 

'’d * [■ j. '’i • . 
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» 
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► 


ITEST = 0 or 1 
Thold = 

Bias » b 


Betad » dB 
ITRNK (normalized) 
MTEST= # of trials 
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Figure 32. Sequential Detection Simulation Flow Chart 
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Figure 34. Variation of Average Synchronization Time and Probability of 
False Alarm With Truncation Time (Normalized) 
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Figure 35. C/A Code Acquisition Time as a Function of Signal-to-Noise 

Ratio (ITRNK= 10) 
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Fii)ut'o 36. C/A Code Acquisition Time as a Function of Siqnal-to-Noise 

Ratio ( I TRNK “ 16) 













b. Variable Detection Points 



Variable 
\ri chip 


NOTE: Position 1 (and hence Signal 1) varies between 0 and 1/2 chip 

randomly. Positions 2,3,4 (and hence Signals 2,3,4) follow 
Position 1 proportionately with constant displacement. 


Figure 37. Two Models of Correlation Peak Detection Used in Simulation 



Table 19. Random Starting Points 
for Monte Carlo Simulation and Resul 




Number of divisions between 
0 and 1/2 chip (N) 

Number of trials (L) 

P(l), which corresponds to 
random position of s(l) at 1 

P(2), which corresponds to 
random position of s(l) at 2 

P(3), which corresponds to 
random position of s(l) at 3 

P{4), which corresponds to 
random position of s(l) at 1 

P(5), which corresponds to 
random position of s(l) at 6 

1 L 

Average P^ = l ^ ’’i 
i=l 
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In Table 20 , the initial condition is picked so that all the 
possible initial values are chosen equally and the average value of Pp 
is calculated, resulting in of 0.9442. The above two cases are com- 
pared to the case of using four fixed signal positions, shown in Figure 
3>a, which has a probability of detection of P^* 0.936. Thus, it is 
(.oncluded that a decision to use a fixed-signal position is preferred 
lOr simulation purposes since the difference between this case and the 
uctual random case is not significant. 

Using the data for the fixed four-signal case, the average syn- 
chronization time, Tg^+T^ 2 » plotted versus SNR, as shown in 

Figure 38, where both synchronization time and Pq are improved over those 
o‘ Figure 35, which corresponds to the single-point case. From Figure 38, 
Pg= 0.983 at SNR = 31.6 dB-Hz. 

Since Pq= 0.9 is a reasonable enough value at SNR =31.6 dB-Hz, 
a refinement of the results is done by picking different values of b and 
so that Pq is within 0.9. The values of b = 3 and 0= -0.9 are chosen 
as optimum values from Figure 39. The corresponding values of probability 
of detection and probability of false alarm are Pq=0.9 and Pp^= 0.001, 
respectively. The corresponding values for as a function of SNR 

is shown in Figure 40. As expected, it is improved by approximately 10% 
at a = 8 dB. 

The bame reasoning is followed at SNR=34.6 Hz, which corresponds 
to ttie C/A code, and the optimum values are taken from Figure 41 to be 
9= -0.07 and b=4.2. The corresponding values for and Pp versus 

SNR are shown in Figure 42. 

6.2.7 Doppler Effect 


The effect of the code doppler on the signal-to-noise ratio a is 
taken in' account in this section. It was shown [4] that the signal-to- 
noise .atio will be changed for each new sample to take the form: 



where otp is the SNR at t=0 (this corresponds to SNR at the peak of the 
triarnlc); Pq is the initial phase error it will take for the four fixed- 
si>-’nal case, 0.75 chips for S(l) or S(4) and 0.25 chips for S(2) or S(3); 


Table 20. Equally Weighted Starting Points 
for Correlation Detection with Resulting Probability 


Number of divisions between 
0 and 1/2 chip (N) 

Number of trials (L) 


P(l), which corresponds to 
random position of s(l) at 1 

P(2), which corresponds to 
random position of s(l) at 2 

P(3), which corresponds to 
random position of s(l) at 3 

P(4), which corresponds to 
random position of s(l) at 4 

P(5), which corresponds to 
random position of s(l) at 5 

L 

i=1 


6 

5 

'’o 

0.944 

0.926 

0.927 

0.959 

0.965 


Average = L 


0.9442 



Average Synchronization Time, sec (T^^+T 




SNR, dB (a) 


Figure 38. Variation of Detection Probability and Synchronization Time 
With Signal-to-Noise Ratio for Four Possible Correlation 

Detection Points 
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Threshold (0) 


Figure 39. Variation of Synchronization Time With Bias b and Threshold 0 for 
Picking b and 0 With Multiple Point Correlation Detection 
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Figure 40. Variation of Pp and With a for b and Optimized 

at .1 = 8 dB for Four-Point Correlation Detection 
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Figure 42. Variation of Total Average Synchronization Time With 
(Also C/N_) for b and B Optimized for a= 11 dB 





and AP is the code doppler, which is assumed in our case to be ±5 chips/sec. 
The effect of the code doppler on the total average synchronization time 
Ts^otal considered for the four fixed-signal level. 

The relation between doppler and and degradation is 

plotted in Figures 43 and 44 for the previously selected values of b 
and 0 at the two levels a=8 dB and a =11 dB, respectively. It can be 
seen that, for doppler offsets greater than approximately 5 chips/sec, 
the performance degrades rapidly. 

6.2.8 Predetection Bandwidth Effect 

The simulation bandwidth was varied among the three values 100 Hz, 
200 Hz, and 400 Hz to determine the performance sensitivity to bandwidth. 
The relation between SNR (dB-Hz) versus was considered for the 

three cases using the best operating points at (peak of triangle) a=8 dB 
and a =11 dB, respectively. The six curves are drawn in Figure 45 with 
the LinCom result, which corresponds to a bandwidth of 500 Hz. It is 
noted from the curves that the Axiomatix results give somewhat better 
T5totai those of LinCom for the C/A code. 

6.2.9 Conclusions 

As a result of the analyses and simulations conducted by Axiomatix, 
it has been shown that a Shuttle GPS receiver with sequential search can 
search and acquire the entire C/A code in approximately 10 seconds. Fur- 
thermore, 1000 chips of P code uncertainty can be searched, with high 
probability of detection, in approximately 16 seconds. An unfavorable 
doppler offset of up to 5 code chips per second has negligible effect on 
the acquisition time. This corresponds to a LOS range rate uncertainty 
of 5000 feet per second for C/A code acquisition and 500 feet per second 
for P code acquisition. 
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Figure 44. Effect of Doppler on Detection Probability and Total Average 

Synchronization Time 
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7.0 Gi^S CARRIER ACQUISITION AND TRACKING 


Most of the GPS user equipment designed to date employs an Auto- 
matic Frequency Control (AFC) loop to provide fast carrier acquisition 
and a Costas loop for carrier tracking during demodulation of the data. 
The composite AFC/Costas loop design provides two major advantages over 
a simple Costas loop: (1) the carrier can be acquired without sweeping, 

and (2) false lock due to data sidebands is avoided. 

To form an AFC discriminator for the loop, it can be seen (ignor- 
ing noise) [5], that an unnormalized discriminator can be formed by 


U) 


I dQ . Q ^ 

^ dt ^ dt ’ 


(73) 


where I and Q are the in-phase and quadrature phase components of the 
signal, respectively. Approximating the time derivatives by finite time 
differences of i, then neglecting the scaling by l/i, (73) becomes 

CO = I(t- i)Q(t) -Kt)Q(t-T) . (74) 

Thus, for a frequency error Aw and a sinusoidal signal, 

I(t) = cos (Acot+4>) 

Q(t) = sin (Ault + 4>) . (75) 

Substituting these values into (73) gives Aco and into (74) gives 
Aui = cos (AvOt - A.oi + >{>) sin (\uit+s’’) - cos (Aot+v,*') sin (Acot - AuU + iji) 

= 2 [sifi (2Au't - Awl + 2i|) + sin (Awi ) - sin (2Awt - Auh + 2^*1 ) - sin (-Aioi )] 

= sin (Awl ) . (76) 

Therefore, using finite time differences, a control proportional to fre- 
quency error can be obtained if |Aau| • 1. Note that (74) is not affected 
by BPSK data modulation unless a data bit transition occurs within the 
delay i. 

If noise is present as in tie GPS receivers, the mean of Aw is 
still given by (76), but there is a SNR penalty due to the multiplication 
of wideband noise processes in (74). To avoid this severe penalty, the 
delay r can be realized by lowpass filters. Figure 46 illustrates a 


composite AFC/Costas loop where the AFC discriminator described by (74) 

Is shown within the dashed lines. 

When the composite loop Is out of lock, a radian frequency error 
Aa>(t) exists In the phase detector output signals £ (t) and c (t). Since, 
typically, the AFC loop bandwidth Is small relative to the noise bandwidth 
of the RC arm filters, then the average AFC output signal <WQ(t)> can be 
computed, assuming Au(t) fixed equal to Aai. In particular, 

S{m(t)m^(t;/W,(j>) sin (Ao)t + (l») 

- m(t) m^(t;Aw,(}») cos (Awt + (|))} , (77) 

2 

where m(t) is the BPSK data modulation, Is the reference signal power. 
Is the In-phase and quadrature phase detector (multiplier) gain, S Is 
the signal power, and 

mj(t;Aa),4)) = 6(p)e^(t) = 6(p) |m(t) sin [Acot + c|)]} 

m^(t;A(jo,(J>) = 6(p)e^(t) = 6(p) {m(t) cos [Aa)t + 4)]} . (78) 

Note that the differential equation resulting from filtering the signal 
by the RC filters 6(s) are written in compact form by introducing the 
Heaviside operator p^d/dt. 

The digital modulation m(t) may be represented as 

00 

'"(t) = I a p(t-nT), (79) 

n=-oo 

where a^ = ±l is the data symbol sequence and p(t) is the data pulse shape, 
e.g. , for NRZ coding, 


(1 for -T/2 < t < T/2 

P(t) = 

10 elsewhere . (80) 


Substituting (79) into (78), 



, .A 










in (t;Aa>,(j») » I a cT~^{G(jw)cjrCp(t- nT) sin (Auit + <|))]} 
n*-® 


J_ 

2j 


I a I a p.Ct - nT) 

■> ** ^ 


n=-oo 


n»-® 


where <?r denotes the Fourier transform with its inverse and 


(81) 


Pl(t) = 


rG[j(u,+ Au))]P(ja,)ej^^^ 


Po(i^) “ [ 6 [j(w - ^)] P(jo)) ^ . 

i-oo 


( 82 ) 


Similarly, 


m^(t;/W),4>) = 2 


I a p^(t-nT)e^^^^^‘**^+ I nl) 

= -OO p = -0O 

( 83 ) 

Substituting ( 79 ), ( 81 ), and ( 83 ) into ( 77 ) and ( 78 ) gives 

n ^ ^ a^a^P2(t-nT)p(t-mT) 


«o(t) = k,\ 2 s 


2 j ^ 

^ m=-® 


CO 00 


"H ^ ^ a^aj^p^(t-nT)p(t-mT) , ( 84 ) 

^ n=-« m=-« 


Taking both the statistical and time averages of ( 84 ) and recalling that, 
for random data. 


a a =6 
n m mn 


1 ; m = n 
0 ; m f n , 


( 85 ) 


_ 1 1 
~T~2~ “ 2i 


Tq/2 


'-Tq /2 n=-« 


I P(t- nT)[p 2 (t- nT) - p^(t- nT)] dt 


we get 


. k -r.-. B 
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and 


<WQ(t)> 

___ 


X lim I [ [ 6[j(«- Au)] P(j<j) P(-ju)' ) 

j f)S«oo ^.*00 

sin [(»-“')V 25 ..j(W)nT(l«dB' 

(-b')V2~' 

-i 11m iff GIj(u + 4»)]P(ju) P(-ju') 

Co ''•CD *»0O 

sin [((.»-(o')Tq /23 j(y.y*)nT dw dm' 


(m-0)’)TQ/2 

Recognizing that 

r ,-j(«-u>')nT . 27T r ./ . 2iik\ 

I e ''' T ^ ^r‘“ ’ 

n=-‘» k=-“> 

then (86) simplifies to 


2ir 2ir * 


(86b) 


(87) 


<WQ(t)> 

~7~T~ 


^ lim I f G[j(m-Au>)3P(jm)Pr.j(u.-^)l 

-■[¥( t )] 


du) 

27T 


¥(^) 

- ~ lim I [ G[j(a>+ Am)] P(jm) p[~-j(w 

J J ^ <30 |(=»oo ^*00 L- — 1 


¥(^) 


do) 

2it 


( 88 ) 


Since the only term in (88) which yields a nonzero contribution is that 
corresponding to k=0, then (88) reduces to 
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<v^(t)> 
'2 2 ■ 


f S^(«.) I« {G[j(«-4,)]) If . 


or, finally. 


where 


. <wITt7> />« 

f(Aui) - — 5 — 5 — ® S (u+Au) Im {G(jw)} 


s>) s 


(89) 


dm 

2ti ’ 


(90) 


(91) 


is the power spectrum of the modulation m(t) and Im{6(j(*»)) denotes the 
imaginary part of the arm filter transfer function. Since G(s) is an 
RC filter with transfer function 


‘ r. j(«/a,^) “ — 


-j(u/w_) 1 

I ’ "c " RC 


(u/u^)' 


then 


Im {G(jo))} = 


oj/o). 


1 + (u»/w )^ 


(92) 


(93) 


Furthermore, for NRZ coded data, 


S_((;o) 

m 


iW 


(94) 


Substituting (93) and (94) in (89), the desired disciminator character- 
istic is given by 



p 


t 


t 


t 


t 


t 


p 


t 




123 


f(Aw) 



— COS AojT + sin AiljT 
“c 



Au) 




Figure 47 illuf *ates f{Aoj) versus Aw/w for w T= 0.837 t (the 
optimum ratio of arm filter 3 dB frequency to data rate for small ST/Nq; 
see [6]) and uj^T=1.66Tr, respectively. The equivalent linear gain, 
of these discriminator characteristics is indicated on Figure 47 by dotted 
lines and is found analytically as follows. Assuming Aw/oj to be small, 

V 

then from (95), 


f(tm) . (96) 

c 

where 



Figure 48 is a plot of versus f T = (d T/27t. Note from Figure 47 
that, as o)^T is increased, the linear gain increases but the region 
of linearity decreases. However, the region of linearity in frequency 
increases as is increased. For example, if the region of linearity 
is defined as an error of less than 10%, then for w^T = 0.83Tr, the region 
is Ad) =0.68 u) or Af= 0.282/T, but is Aw = 0.446 w or Af= 0.370/T for 

W \0 

w^T = 1.667T. Figure 49 plots the region of linearity versus f^T. 

Since is approximately 1 for f^T>.10, the region of linearity is 
approximately 












7.1 AFC Loop Response - Linear Model 

We observe from Figure 47 that, for a given value of <u^T, vhe 
d'iscrimlnator characteristic Is linear over a certain range of frequency 
offset and, In particular. Is described by (96) In this recion. Moreover, 
since (95) and hence (96) were derived under the assumption of a constant 
frequency offset, I.e., o)(t) - w(t) = then from Figure 46, the closed 
loop AFC response is described by 


«i(t) = KyF2(p) K^\^SK^(oi(t) - <i(t))^u,^ (99) 

or, in Laplace Transform notation. 



SKK^ 

1 + SKK^F2(s)/u.^ 


(100) 


A 2 2 ★ 

where Ky is the VCO gain in rad/sec/v and K-KyK^^K^^. Assuming that 
F 2 (s) is an imperfect integrating filter, i.e.. 


then 


F2(s) 


1 + ST 
1 + S T 


22 

12 


( 101 ) 



S K (1 + s ^22^^ ***c 
1 + S K Kfjp/ <0^ + s(t-j2‘*’SK K^t 22 / ) 


( 102 ) 


which implies a first-order loop response with radian cutoff frequency 

1 S K K^/u>^ 


(tt 


nF 


t^2 ^ SKK^T22/<*>^ 


(103) 


K(h is not included in the definition of K since it is actually 

» ^1 - A . ^ 


a function of ai_T rather than a constant. 
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If a conventional Costas loop were to employ F 2 (s) of (101) as its 
loop filter, then the loop response would be second-order and character- 
ized by the following parameters: 


rg ^ SKD|^t22F()2 ' 



^22 

^2 


“n2 



(104) 


In (104) w ^2 the natural frequency of the loop, and for small Fq 2 (the 
usual case in design), 

rj i . (105) 


where ^2 the loop's damping factor . Also in (104), 


0 


m 


4 




( 106 ) 


is the data modulation distortion factor (equivalent power reduction due 
to arm filters).* Using these definitions in (102), we obtain 



Ihe function ) which appears in many places in (107) and (108) can 

be evaluated as follows. For an RC arm filter and NRZ data, one can show 
from (106) that 


* 

The parameter 0^^^, like , is a function of (jj^T. 
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009) 


Thus, from (97) and (109), 
K.v / K 


with 




. y V 0)J-2+(2 + u)J) 6“*"*=^ 


w^T - 1 + e 
c 


■Uj-T 


( 110 ) 


( 111 ) 


For example, at the small ST/Nq optlntum design ratio of arm filter eutoff 

frequeney to data rate, viz., fa»^T*0.83TT, then K^* 0.216. Figure 50 

plots K. versus f_T. Note that, for ri>10, K «l/o) T. 
c c c c c 


7.2 AFC Loop Response to an Input Frequeney Step 
Using (110), we ean write (107) in the form 


u)tS 

U) ' s1 


'‘cT 


r “1 


ri BC -1 

A + Bs 

s T 

* • 6 

C + Ds 


Lc-Os “J 


where 


A = (1), 


n2 


■fl", 


h2 ’ 


(112) 


C ■ ^ ‘ “ * ’ * "c ^ “n2 • 

■n 

If w(t) Is now a frequeney step AQu(t), i.e., w(s) * Aq/s, then 


<i(s) » K^l 


^0 - f\ 

^ \C + Ds/ 


^ 1- 

'' ' + 0 0 


^e^T^O 


1 B/D 
' ' A/C 


0 . 


*c/»J 


(113) 


(114) 
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From (113), 

1 1/D 

‘ ■ A/C 



1 


1 



(115) 


since simplify (114) any further, we must express 

in the denominator of (115) in terms of the parameters of the conventional 
Costas loop portion of the overall composite loop. Assuming for F^(s) an 
imperfect integration filter of the form 


F^(s) 


1 + S T 
1 + S T 


11 


(116) 


then the Costas loop is second-order with parameters [analogous to (104)] 


1. 


r, ^ SKD„To, F 


2 ] 
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m ‘21 ' 01 ’ 


01 


^11 




U) 


hi 


'21 


'Ll 


VI 
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(117) 


where, in addition, is the single-sided bandwidth of the Costas loop 
portion alone. Thus, from (104) and (117), we have that 

•^c ^ 'J'^2 “'n2 " *^0 ^ ^ *^m *"o2 


'21 


/'jD2\ 

t 

(l02\ 

'^01 ' 




Since a conventional Costas loop is typically designed tor critical damping 
(r^ =2) and B|^^T «.< 1, then if we arbitrarily choose Fq^ ^ F^^ , we have 
that 


^c^ ‘ ''2 '"n2 



(119) 


or equivalently, from (113) and (116), 

nil- 1 - _D/p_ ; 

^ ' • ' A/C 


(120) 



132 


< 1)1 SJ c 


(121) 




'c ' 2 “’n2 


where we have defined 


( 122 ) 


8Kc T 


(123) 


The equivalent time domain response is then 


«(t) - [t - u(t) . 


(124) 


8 T 


(125) 


The corresponding frequency error response, = w(t) - o)(t), is given 
by 




(126) 


Furthermore, since from (121), the AFC acts like a first-order loop, its 

two-sided noise bandwidth W, is related to <o ^ by 

LF nF 

. , „ 2 ;,) 

Lr 1 ‘ I ^ 
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Alternately, the single-sided noise bandwidth Bj^p = 's related to 

n _ “nF 


“nF 


LF 


(128) 


Substituting (128) Into (126), we obtain the desired radian frequency 
error response: 




(129) 


7.3 Acquisition Behavior of the Composite AFC/Costas 
Loop In the Absence of Noise 

It Is well known that an Imperfect second-order phase-locked loop 
is characterized by a finite acquisition range (often called capture or 
pull-in range) which represents the largest frequency difference between 
the frequency of the received carrier and that of the VCO for which phase- 
lock will occur. Typically for high-gain loops, this acquisition range 
(relative to the loop-bandwidth) is solely a function of the loop damping 
and the ratio of the time constants of the imperfect loop filter. In view 
of the similarity between the differential equations of operation of the 
conventional Costas loop and the phase-locked loop, one might expect the 
two to have similar acquisition behavior. Since, in effect, a Costas loop 
reconstructs and tracks a phantom carrier at twice the input frequency, 
the above statement is indeed true, at least for frequency errors on the 
order of the loop bandwidth. The similarity to a phase-locked loop, how- 
ever, tends to break down as the frequency error becomes large, due to 
the tendency of the Costas loop to false lock to a sideband of the data 
modulation. Such a false lock phenomenon does not occur in a phase- 
locked loop, since it tracks an unmodulated carrier. This false lock 
tendency of the Costas loop often prohibits acquisition by sweeping the 
VCO frequency over the range of frequency uncertainty, particularly if 
the false lock frequencies, which occur at integer multiples of half the 
data rate, fall within the uncertainty region. 
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The addition of the AFC network to a conventional Costas loop 
considerably improves the false lock problem and increases the acquisi- 
tion range to a value on the order of the arm filter bandwidth. 

To begin the analysis of the acquisition behavior of the composite 
AFC/Costas loopt we present (by analogy to the phase-locked loop) the 
results for acquisition range and frequency acquisition time* for the 
Costas loop portion alone. 

For a phase-loc''*‘d loop with frequency error Af, the loop will 
lock if [7]: 


l-'fl ^ (’30) 

and the corresponding time to achieve the frequency lock is [7]: 



Since the differential equation of operation of the Costas loop is iden- 
tical in form to that of the phase-locked loop except that is replaced 
by 2(J) and Af by 2Af, then (130) and (131) apply for the Costas loop if 
Af is replaced by 2Af. In particular, for the Costas loop portion of 
the composite loop of Figure 46, we have that frequency lock is achieved 
if 


|Afl 



and the time to reach this condition is 


(132) 


We shall ignore the phase stabilization time (the additional time 
to reach phase lock once frequency lock has been established) since this 
time has only a small effect on acquisition time for appreciable initial 
frequency errors. 
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(133) 


where r^, and are the loop parameters defined in (122). 

Since the AFC increases the acquisition range by producing a 
restoring force proportional to the frequency error [see (95) or its 
linear equivalent (96)], it is essential that the value of initial fre- 
quency error be such that the AFC can reduce this error to a value within 
the pull-in range of the Costas loop. Since, from (129). the initial 
radian frequency error can, in the limit as t approaches infinity, be 
reduced no further than Aq 5/(1 + €), then for acquisition of the composite 
loop, we require that 



(134) 


Letting r^ = 2 and using the corresponding definition of C given in (125), 
we arrive at the inequality 



(135) 


Since, from Figure 47, the AFC discriminator characteristic is 
approximately linear over the frequency range AuxO.68 cj for <d T = 0.837t 
with slope then the maximum acquisition range of the composite loop 
is approximately 


lAf| 1 0.282 

BlI 2. • 


(136) 


If Aq is set equal to 0.68 and F^^ = 0.001, then at = 0.83tt, the 
inequality in (135) can be reduced to 

> 0.00637 . (137) 


I 



The corresponding maximum acquisition range for the Costas loop alone 
and the same loop parameters Is obtained frcxn (132) and Is given by 


/iM\ 

), 


= 9.488 


Ll /max 


( 139 ) 


Thus, a modest Improvement In acquisition range is possible through the 
addition of the AFC. More Important, however, this increased acquisi- 
tion range Is accompanied by a considerably reduced acquisition time. 

When the frequency error is large, the behavior of the composite 
loop Is controlled primarily by the AFC loop. Thus, the time t^ to 
reduce the frequency error from Its initial value AfQ-AQ/27r to another 
value, say Af^, can be determined from (128), viz.. 


2’TAfl = 2irAfp(^^J + 

2TTAfo 

1*5 ® 

\140) 

^1 “ ‘ 

r^f, 


(141) 


When Af^ becomes within the pull-in range of the Costas loop, then the 
time t 2 to complete the frequency acquisition process is determined by 
the combined action of the AFC and Costas loops. This time will, in 
general, be less than that computed for the Costas loop alone attempting 
to achieve frequency lock from an initial frequency error of Af^. The 
latter is computed from (133) with Af replaced by Af^. Since an expression 
for t 2 is difficult to obtain, we shall upper bound t 2 by the value 
computed for the Costas loop acting alone as discussed above. Thus, we 
propose to upper bound the acquisition time T^^q of the composite loop 
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by the sum of as determined by (141) and t^ as determined by (133) with 
Af^ replacing Af. Since Af^ is thus far arbitrary, we shall choose its 
value so as to minimize = tj + tg . Furthermore, if Af^ is consider- 
ably less than the pull-in range of the Costas loop, then the acquisition 
time formula of (133) simplifies to 


ri.l 


4 B 


LI 


/r, + K /2Af, 




(142) 


We shall make this assumption in what follows and check its region of 
validity later on. Since we wish to choose Af, so as to minimize 

I * 

Tacq “ ■*■^2 ’ differentiate Taqq with respect to Af^ and equate 

the resulting expression to zero. Letting 


2 I r, / 


C = 


and 


X = 2Af, 


we have from (141) and (142) that 


(143) 


'ACQ 


4B 


In 


LF 




2a f. 




(144) 


Thus, 


9T 


M = 

3 X 


1 


4 BlF X - ^^^0 ^ 


1 + C 


or, after simplification. 


3 ^^^0 ^ 2 

X - , . - X - 


1 +? 


4C‘ 


(>"i ■^1) r 2 in 

* Arc" * ‘ ° 


2CMr^ + l) g 


LF 


Ll-l 


X + 


AfflC 


2C^ (1 + 4) 


(145) 


= 0 . (146) 


3 2 

A cubic equation of the form x + px +qx + r = 0 may be reduced to 


y + ay + b = 0 


( 147 ) 


by substituting for x the value of y-p/3. Here, 





a = q 




b = * 2 ^ (2p^ - 9pq + 27r) . 
The case of interest here is where 


(148) 


h2 ,3 

-g- + ^ < 0 , 


(149) 


for which there will be three real and unequal roots. These solutions for 
y are obtained as follows. Compute the value of the angle e in the 
expression 

cos 0 = . (150) 

r~jT 


(150) 


Then the three solutions y^ , y 2 » and y^ are 


y, ■ 2 


a 6 

T 


^2 = 


' 3 

COS 

( t * t ) 


(151) 



1--^ 
' 3 

cos 

( t ^ t ) 



equivalently. 

the 

three 

solutions for 

X are 


^1 


/-f 

COS 

± 

3 3 



X2 = 


l-f 

cos 

( t * t ) 

. -B- 
3 

(152) 
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l-f 

cos 

(f^r) 

. ^ 

3 • 



Comparing (146) and (147) and normalizing x by 2AfQ, we find that, cor- 
responding to (152), we have the three solutions 
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Although ft appears from (153) that three solutions exist, In reality 
only the first, namely ' -j/ZAf^, results In 


0 + 0 



0 . 


( 157 ) 


The Inequality of (157) Is necessary to Insure a real solution for T^j.q, 
since we observe that the above expression Is the argument of the natural 
logarithm In the first term of (144). Finally, then, solving the first 
equation In (153) using (154) and (156), we can compute our desired upper 
bound on from [see (144)]: 



As an example of the acquisition time for a GPS receiver, consider 
the parameters used In the GPSPAC. The data rate Is 50 bps and = 35 Hz 
(or Bj^^T = 0.7), '"oi .001, and r.j - ^ For optimum tracking at low ST/Ng, 

Figure 51 illustrates the acqul- 

From 

Figure 51, it Is seen that there is little difference In the curves for 
the two values of f^T, but the region of linearity as indicated by the 


the arm filter bandwidth is f^T= 0.41b 

sitlon time with these parameters with fT= 0.415 and f T*10. 

c c 


dashed lines of 


is larger for f^T = 10. 


Also Included In Figure 
51 is the acquisition time performance of a conventional Costas loop as 
computed from (133) which assumes » 10. The asymptote Indicated 

by Afg/Bj^^ = 9.488 is the pull-in range of the Costas loop determined from 
(132) at which point the denominator of (133) vanishes and the acquisition 
time becomes infinite. The acquisition time for the Costas loop and the 
composite AFC/Costas loop is ^^(;q" 13.7 sec and T'j^cQ"^3.4 msec, respec- 


160 Hz, f^ = 500 Hz, and Bz. 


In other 


tively, for 

words, the composite AFC/Costas loop acquires 316 times faster than the 
conventional Costas loop. These acquisition times also assume that the 
conventional Costas loop does not false lock; but, there are false lock 
points at ±25, ±50, ±100, and ±150 Hz about the carrier frequency for the 
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conventional Costas loop and no false lock points for the composite AFC/ 
Costas loop. 

Narrowing the Costas loop bandwidth Improves tracking under low 
dynamic conditions such as the Shuttle/GPS receiver will encounter on- 
orbit. Figure 52 presents the acquisition time for the AFC/Costas loop 
with B|^iT»0.1 or seven times narrower than In Figure 51. For Hz, 

the acquisition time of the composite AFC/Costas loop Increases by a factor 

of 430 with ® ^ACQ“ " \f* 

In fact. Figure 52 shows that, for f T» 10, the acquisition time Increases 
exponentially as Afg Increases above 3 Hz for ^^ 2 . However, 

the composite AFC/Costas loop can still acquire at 160 Hz, while the con- 
ventional Costas loop can never acquire. If narrower arm filter band- 
widths than fj,T*10 are used, then as shown In Figure 52 the acquisition 
time Is decreased within the region of linearity for the narrower band- 
width; but, of course, the acquisition range Is decreased. 

One final point before leaving the discussion of acquisition per- 
formance regards verification of the assumption made just prior to (142), 
allowing use of this simplified acquisition time formula rather than (133) 

In deriving the bound on acquisition time. The assumption made was that 
the solution found for the value of Af^ which minimized would satisfy 



or, equivalently. 



(160) 


Usinq x^/(2Afp) as determined from (153), Figure 53 plots the 
left-hand side of (160) versus for r^ = 2 and Fq^* 0.001. For 

these parameter values, the right-hand side of (160) equals 9.488. Since 
the denominator of (133) Is computed as one minus the square of the ratio 
of the left-hand to right-hand side of (160), we observe from Figure 53 
that this correction factor is well within our approximation over most 
values of frequency offset within the acquisition range. 






Figure 53. Check on Validity of Using Simplified Formula for Calculation of Costas Loop Acquisition Time 
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Tracking Performance of Composite AFC/Costas Loop 


When separate loop filters are employed for the AFC and Costas 
loops as in Figure 46, the equivalent additive Costas loop noise and AFC 
error signal (which acts like additive noise when the loop is locked) are 
processed by different closed loop transfer functions insofar as evaluating 
phase tracking jitter. To see this, we begin by relating the VCO output 
phase to its two input signal components, i.e.. 


§(t) = j [F,(p)z„(t) ♦ F2(p)Wj{t)]. 


where ZQ(t) and WQ(t) are given by 


z^it) = {[v^m(t) - s,n 2^(t) 


+ - Ns(t)3 cos 2*(t) 


(161) 


(162) 


iNc(t) [*^m(t) - N^(t)] - f)^(t) [v^m(t) - N5(t)3} 

(163) 

where the additive channel noise is expressed in the form of a narrowband 
process about the actual frequency oi the observed data, viz., 

n^(t) = *'^|N^(t) cos [(Ogt + 0(t)] - N^(t) sin [wQt + G(t)]|. (164) 

The narrowband processes N^(t) and N^(t) are statistically independent, 

stationary, white Gaussian noise processes with single-sided noise spectral 

\ \ 

density Nq W/Hz. The narrowband processes N^(t) and N^(t) in (162) and 
(163) are N (t) and N (t) passed through the lowpass arm filter 6(s), 
respectively. Alternately separating ZQ(t) into its signal and noise 
components, we can rewrite (161) as 


§(t) = JL_^ p^(p) + -J- V2(t.2*)J + F2 (p) w(t) . (165) 
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V 2 (t, 2 ^) ^ t-N(f(t) + n^(t) - 2 /r m(t)Ng(t)] sin 24 . 

+ [ 2 v^‘ m(t)N^(t) - 2 N^(t)N^(t)] cos 2 * (166) 

and 

w{t) = 2wQ(t) = 2^m(t)N^(t) - 2N^(t)Nj(t) 

- 2v^ m(t)N^(t) - 2N^(t) Ng(t) . (167) 

Then, letting K - replacing m^lt) by - <m^(t)>, and ignoring 

the self-noise of the modulation, we get 


2e(t) = y {F^(p) [S D|^ sin 2((.(t) + V2(t,2^)] + F 2 (p) w(t)} . (168) 

Finally, letting <j> 2 (t) = 24 i(t) = 2[e(t) - e(t)] = ©^(t) - § 2 (t) and taking 
the linearized case, i.e., sin 2if(t) = 24i(t), gives the desired differential 
equation of loop operation, viz., 

<^2(t) = 02(t) - y F^(P)S0^ ^ ''2^^*‘^2^ ^ 


or 


4’2(^) 


p + Sd|^KF,(p) “2*'^' * 

SP„KF,(p) 

P + SD^KF,(p) 




(P) 


sDjrri^ 




S D. 


m 


)( t ,4>2) 

Td~ 


(170) 


Defining, as is done for conventional Costas loops, the closed loop 
transfer function. 


SD KF|(p) 

”l^P^ " P + SD^KF“(p) ’ 


(171) 
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equation {170) simplifies to 


Vo(t»4p) 

i>2(t) = [1 - H^(p)] 02(t) + + 

m 


k(p) 


F,(P)' 


TOT 




SD. 


m 


= [1 - H(p)3 02(t) + Ng(t) 


(172) 


where N (t) is the total equivalent noise and, as predicted, is composed 

V 

of the sum of the additive Costas loop noise V2(t,24>) and AFC error signal 
w(t), each processed by different closed loop transfer functions. To 
evaluate the tracking performance of the composite loop, we must find the 
spectral density of Ng(t). 

From the definition of N^(t) given in (172), we see by inspection 

that its power spectral density Sm (w) is given by 

"e 
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H^(jw): 


|Fp(ju.)|2 



7 Fp(jo)) 

n — S (ui) 
F^(jo)) Vpw' ' 



F^(ju.y ^WV2^“^ 


(173) 


where the asterisk denotes complex conjugate, Sy 2 (“) and S^(w) are the 
power spectral densities of V 2 (t,({i 2 ) and w(t), respectively, and Sy^^(o)) 
and S^y^(u) are the corresponding cross-spectral densities. Since the 
bandwidth of the Costas loop is ordinarily designed to be narrow with 
respect to the equivalent noise bandwidth of v. 2 (t,({> 2 ) and likewise for the 
AFC loop with respect to w(t), we can approximate V 2 (t,i|i 2 ) and w(t) as 
delta-correlated processes with power spectral densities Sy^(O) and S^(0). 
and cross-spectral densities Sy^w(O) and From [8], 

Similarly, from the definition of w(t) in (167), we find after some alge- 
braic manipulation that 
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SJO) = 16s(^)| S^(o)) [Im {G(ja,)}] 

00 

+ 16(^)^j [Im {G(ja,)}]^ 


2 ^ 
2tt 


V°' “ ' “ 


(175) 


Thus, substituting (174) and (175) into (173), we get 


e (S D^) 


+ ' — ^ 




Where 

K. * 


. rje(>)l^ I; 
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K, 


F2(jo)) 1 2 


F, (ju)) 


t2SNoD„,A'p + NQ^B.^J, (176) 




r.s^(^) [i» iG(j.o)i)^ ^ 
r.[Im (G(>))]^f 


/.:Nj»)i^ if 
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s„(-) |G(jo.)|^ g ; 8, = J |G(j»)|‘^ 


2 du) 


(177) 


The total equivalent noise power, o|^ , can now be obtained by integrating 
(176) between and ®, viz.. 


2 L 


s„ (u.) » 

\ (SD^j 


rrT[2SN„D„K„*N„2B,K ]f" 

N n ) ^•oo 


,00 


(SD„) 
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. (178) 
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The first integral In (178) is the well-known relation for loop bandwidth 
of a second-order Costa** loop* i.e.. 


r 




(179) 


The second integral in (178) is evaluated as follows. Define 
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F2(j<‘>) 

iHi(ju.)r 

4m 00 

f^TW 
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^00 


where 


F«(s) 


A^(s) - H^(s) 


(180) 


SDn,K.s[SD^K(T33.T^^)].s-^SD^KT^^t3g 

S Dn, K+ s[S Djjj K(t2^ + t^2) + j + t^2(S K T2^ + 1 )] + s^ 


^ * ^2^ 
djj + d^s + d2S^ + 


(181) 


Integrals of the type given in (180) have been previously tabulated (see 
for example, [7], p. 135, Table 4-1). when the denominator of Aj^(s) 
is a third-order polynomial, i.e., n = 3, then the result is 
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3 




* (cf -gCpCgjdoda ♦ 


(182) 


Substituting (181) into (182), one obtains after considerable simplification 



(183) 


where we have made only the assumption that Fq^,Fq 2 « !• Further relating 
®L12 bandwidth of the Costas loop alone by 
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®L12 “ ^B®L1 ’ 


(184) 


one can show, using (104), (108), (111), (117), (128), and (183), that 
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^c\lT^ -^r— 


(185) 


Thus, for (185) reduces to 


B. .\2 




Bli t 




(1 + c) 


(186) 


where ^ is defined in (125). Figure 54 illustrates the bandwidth ratio Yg 
as computed from (186) versus 8^p/B|^^ for Fq^ = 0.001, Bj^^T=0.7, 0.1, and 
fj= 0.415, 10. 

Finally, combining (178), (179), (183), and (184), the total 

equivalent noise power, which is identical to the mean-square phase 
2 

tracking jitter o<{, 2 , is given by 

_2 _ 4^4 


P S. 


ps, 


2^2 

o, + 0 , , 


(187) 


where 


"o\l 


is the signal-to-noise ratio of an equivalent linear loop, e.g., a phase- 
locked loop, having a bandwidth equal to that of the Costas loop, is 



Figure 54. Bandwidth Ratio Yg Versus the Ratio of AFC to Costas Loop Bandwidth 
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the familiar squaring loss of the Costas loop, namely, 


o _ _ III 

K[- 

* kC X ^ 

“ ‘■i 


and is the squaring loss of the AFC loop which is defined by 


\ 1 

+ K- 

L“ "tM 


In (188) and (189), p^-2S/NqB^ is the signal-to-no1se ratio in the arm 
filter bandwidth. For a single-pole RC arm filter as in Figure 46, 


\ “ T * 


Furthermore, for NRZ data, Kq and Jig defined in (177) can be evaluated as 


3- (3 + co T) e ^ 


1 . 1 [1 . 

w T 

c 


1 "U)^T 

"d = i • (’’■ 

c 

Figure 55 illustrates Kp, and as defined in (109), (191) and (192) 
versus f^T = a)^T/2TT. At f^T = 0.415, (191) and (192) are evaluated as 
Kjj= 0.782 and A'p = 0.437. Also, from (109), 0^= 0.645. Similarly, at 
f^T*10, then Kp = 0.992, Ap= 0.0162, and 0^^^= 0.984. Thus, the component 
of phase tracking jitter due to the Costas loop for f^T= 0.415 is given 


R [o.782 + 
4NoBli 


0.645 


(193) 
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and for f T* 10, 

V# 


2 . 


0.992 + 


O.SOffl 




0.984 


Since for a single-pole RC filter, B^»oj^/ 2, we can express as 






Thus, (193) can be rewritten as 
. 2 ‘""0“U 


“I — r 

and (194) can be rewritten as 


1.2,2 . ,.567(!i»^)(j^j 


,2 . '•VU 

'l - “1“ 


,008 . ,6.223(!!»^)(5j^j 


Similarly, the AFC loop component of phase tracking jitter for f^T 
is given by 



I 


(194) 


(195) 

(196) 

(197) 

= 0.415 

(198) 


2 ( 

(199) 
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Mgyres 56 and 57 Illustrate the mis phase tracking jitter a^2 

9 as computed from (196)- (199) and (187) versus ST/Nq (Ejj/Nq) In dB for 

0.001, r^ = 2, and B|^^T=0.7. 0.1, respectively. The parameter In 
both figures Is the ratio of AFC to Costas loop bandwidth for f^T» 0.415 
and f^T«10. Acquisition may take place If Ojj^gll.O radian. The tracking 

• requirement for the Shuttle/6PS Is 0*524 radian. In Figure 56, the 
composite AFC/Costas loop tracking performance Is the same as the conven- 
tional Costas loop when B^^p/Bj^^lO.Ol for f^T= 0.415 and when B|^p/Bj^^ < 0.001 

for f T=10. In Figure 47, the composite AFC/Costas loop tracking per- 
\0 

* formance Is the same as the conventional Costas loop when Bj^p/B|^^ <.0.001 
for f^T=0.415 and when B|^p/B|^^ £ 10"^ for f^T»10. Comparing Figures 56 
and 57 with Figures 51 and 52, one Immediately observes the trade-off 
between acquisition and tracking performance of the composite loop as the 

♦ Bj^p/Bii bandwidth ratio Is varied. Alternately, after acquisition, the 
AFC discriminator input to the VCO can be switched off and the Costas loop 
arm filter and loop filter bandwidths can be readjusted for tracking. In 
this case, with a little more complexity, the acquisition and tracking 

♦ performance can be optimized separately. In fact, many GPS user sets do 
switch the AFC loop off and change the filter bandwidths for tracking. 
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Figure 57. Tracking Jitter Performance for Composite AFC/Costas i.oop 

(Bj_,T= 0.1) 



8.0 


GPS RECEIVER CLOCK MODEL 


8.1 Introduction 

Oscillators used for timing in the GPS user set and satellite 
are characterized by long-term instability and short-term ins*:abil ity . 
Long-term instability is based on measurement intervals of hours or days 
or longer; short-term instability is based on measurement intervals of 
minutes or seconds or shorter. Short-term instability accounts for the 
random and relatively rapid frequency (or phase) fluctuations that occur. 
Long-term instability accounts for the constant frequency error and 
relatively slow frequency drifts that occur due to aging. This section 
provides the quantitative characterization of oscillator stability that 
is necessary for specifying the Shuttle GPS receiver oscillator. 


8 . 2 Technical Discussion 

A typical sample function of the radian frequency process of an 
oscillator is illustrated in Figure 58. For a square clock signal, the 
oscillator outputs s(t,<f(t)| can he represented as 

s^t,4’(t)^ = Sin 4>(t) 

where 

4»(t) = (i)q 5^t + lii(t) 


and (Oq 
ut 
i(t) 


frequency if the oscillator were perfect 
drift or long-term component 
short-term component or phase noise. 


Thus, 


■t(t) 


i(x) dA 


= WQ(t - tp) + s^(t - Tq)^ + lii(t) - >J'(tQ) 

characterizes the instantaneous phase accumulations. Figure 59 shows a 
typical aging curve for a fairly good quartz crystal oscillator. Least 
squares fits of the straight lines to the two parts of Figure yield agi 
rates of 0.536 ' per day and 0.515 ^ lO"^*^ per day, respectively. 
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For the GPS clocks, both settability error and long-term drift 
t contribute to the long-term elapsed time error, Settability is defined 

as the precision to which two clock standards can be set to the same fre- 
quency. Just how accurately two clocks can be set to the same initial 
frequency is determined by F, the fractional frequency offset coefficient, 
t The accumulated time error due to settability is given by 

aT^ = 8.64 X 10^*^ FT. 

-7 -13 

^ where aT^ is in microseconds, F is less than 10 but greater than 10 

in parts per 10~^ per day, and T is the elapsed time in days. 

Long-term drift arises from slow changes in the average frequency 
of the clock due to changes in the clock resonator. The accumulated 
^ time error due to long-term instability is 

aT^ = 4.32 X 10^^ F'T^, 

where AT, is in microseconds, F' = n/2iT and is less than 10"^ but 

«- _'I3 _v 

♦ greater than 10 in parts per 10 per day, and T is the elapsed time 

in days. The total accumulated time error after T days is 

at = 8.64 X 10^° FT + 4.32 x 10^° F'T". 

f Table 21 summarizes the characteristics of several oscillators 

that are used as clocks. The frequency-time standards are usually of 
1 MHz, 5 MHz or 10 MHz. The oscillator used by the GPS satellites is 
the Efratom Rubidium 10 MHz with characteristics shown in Table 
t The short-term instability is of more interest in the use of the 

GPS than long-term instability because the operation of calculating 
pseudo range and other processes of time in the navigation filter are on 
the order of 0.1 second to a few seconds. Therefore, ifi(t) must also be 
> considered. In general, the radian frequency is characterized by 

Ut) = a + ^(t) 

with a an independent random variable and ^(t) a zero mean stationary 
« random process. The phase accumulated in x seconds is 

ij;{t+T) - i|)(t) = at + i)(t+T} - <{>(t) 


I 





which has a mean square value cf 


D^(t) = 7 + D^(t) 

where the overbar denotes the statistical expectation. Let the average 
phase Instability (or change In phase over a measurement of t seconds) 
be denoted by then 




“0^ 


ojqT 


( 200 ) 


There are several models for the D^(t) of an oscillator which character 

<p 

Ize (1) RC-type frequency Instabilities, (2) Gausslan-type frequency 
Instabilities, (3) flicker-type Instabilities, (4) perturbation noise, 
and (5) oscillator spur and vibration noise. For example, with RC-type 
frequency instabilities 

^ " 1 ]} . 


where a» is the rnis frequency deviation of the oscillator and t,, is the 

<p u 

correlation time. For short-term instability, D^(t) has an asymptote of 
Therefore, from (200), 

/tT 2 

Mil - nlUL!L 

(OqT -OqT 


and the rms elapsed time error due to short-term instability is 



10^ ft , 


where AT, is in nanoseconds, t is in seconds, and f= 1S 

V 9 

given in Table 21 under short-term instability for 1 second averaging 
time. Note for this model of oscillator instability, the short-term 
instabili / for 0.1 second averaging time in Table 21 can be used when 
T is less than 0.3 second. 
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9.0 SHUTTLE GPS INTERFACE CONTROL DOCUMENT 

Because of the criticality of the Shuttle navigation function 
and because GPS navigation entails elements external to the Space Shuttle 
Transportion System (STS), it is necessary to establish an Interface Con- 
trol Document (ICD) between the GPS system and the Shuttle. Examples of 
the external elements include the NDS Space Vehicle Navigation Subsystem 
and the PRN Navigation Assembly signal characteristics. 

Such an ICD has been established for the conventional user sys- 
tem segment, which does not include the Space Shuttle or any other 
space-borne users of GPS. This ICD, MH08-00002-400, Space Vehicle Navi- 
gation Subsystem and NTS PRN Navigation Assembly/User System Segment and 
Monitor Station, is a valid starting point for the GPS/Shuttle ICD. It 
alone, however, is not sufficient to guarantee the succcessful utiliza- 
tion of GPS by the Shuttle. An example of one shortcoming is that this 
ICD specifies the GPS signal power level at the earth's surface but not 
at the orbital altitudes at which the Shuttle operates. The current ICD 
for earth-bound users reflects the fact that the GPS satellite antenna 
patterns are shaped to provide a maximum power density approximately uni- 
formly over the earth's surface, with rapid fall-off beyond the edge of 
the earth. The Shuttle/GPS ICD must specify the equivalent "pointing 
loss" for several Shuttle orbital altitudes or must specify the Effective 
Isotropic Radiated Power (EIRP) versus angle of the satell ite/earth 
centerl ine. 

A preliminary outline of the Shuttle/GPS ICD is given in 
Table 22. This outline will be expanded and the ICD written according 
to the expanded version during the forthcoming follow-on contract. 
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Table 22. Preliminary Outline - Shuttle/6PS ICO 

t 

1.0 

SCOPE 


2.0 

APPLICABLE DOCUMENTS 


3.0 

REQUIREMENTS 


3.1 

Interface Definition 


3.1.1 

GPS Satellite Navigation Subsystem 


3.1.2 

GPS Signal RF Characteristics 

f 

3.1.3 

GPS Signal Baseband Characteristics 


3.1.4 

GPS Time 


3.1.5 

Space Shuttle Navigation Subsystem 


3.2 

Interface Identification 


3.2.1 

GPS System 


3.2.2 

Space Shuttle System 

t 

3.3 

GPS Satellite Constellation Definition 


3.3.1 

Orbital Characteristics 


3.3.2 

Deployment Schedule 


3.3.3 

Maintenance/Replacement Schedule 


3.4 

Radio Frequency Characteristics 


3.4.1 

RF Signal Structure 


3.4.2 

Navigation Data Update 


3.4.3 

GPS Satellite Operating Modes 


3.4.4 

Frequency Plan 


3.4.5 

Out-of-Band Emissions 


3.4.6 

In-Band Spurious Emissions 


3.4.7 

Signal Coherence 


3.4.8 

Carrier Phase Noise 


3.5 

Code Waveform Characteristics 


3.5.1 

Pulse-to-Pulse Jitter 


3.5.2 

Incidental AM 


3.6 

Group Delay Variations 


3.7 

Timing Accuracy 

» 

3.8 

Antenna Polarization 


3.9 

Shuttle Received RF Signal Levels 


3.9.1 

On-Orbit 


3. 9. 1.1 

LI Navigation Signal 


3. 9. 1.2 

L2 Navigation Signal 


3.9.2 

In-Atmosphere 


3. 9. 2.1 

LI Navigation Signal 


3. 9. 2. 2 

L2 Navigation Signal 


3.9.3 

Pointing Loss Slope 


3.10 

Navigation Signal Structure 


3.10.1 

P Code Generation 


3.10.2 

C/A Code Generation 

» 

3.10.3 

P Code Acquisition From C/A Code 


3.11 

Frequency Standard 


3.11.1 

Accuracy 


3.11.2 

Drift 


3.12 

Signal Data 


3.12.1 

Data Waveform Characteristics 


3.12.2 

Data Format Characteristics 


3.12.3 

Data Block Contents 


3.13 

GPS Denial of Accuracy 


3.13.1 

Algorithm 


3.13.2 

Shuttle Navigation Accuracy 

■? ■ • 
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APPENDIX A 


MEASUREMENT OF TACAN TRANSMITTER PULSE SPECTRUM 


The following is a description of the measurement performed by 
JSC of the TACAN pulse spectrum. 

The TACAN used was a Hoffman AN/ARN84 operating on channel 26 
(X mode) in the search mode (150 pulse pairs per second). The test 
equipment used is shown in Figures A-1 and A-2. 

The transmit twin pulses were recorded as shown in Photograph A-1 
using the test equipment configuration as shown in Figure A-1. The pulse 
jitter technique used for coding caused some synchronization problems but 
the pulses obtained are typical of those described in the TACAN literature. 

Photograph A-2, the transmit spectrum, was obtained as shown in 
Figure A-2. The total of 63 dB attenuation was required in order to pro- 
tect the spectrum analyzer from burnout. The photograph shows the main 
spectral component at 1150 MHz and no discernible component at the 
1226 MHz region. Photograph A-3 is the same configuration as Photograph 
A-2 except that the baseline clipper control on the analyzer was turned 
to minimum. This caused the trace to bloom on the scope. The scope was 
readjusted to remove the bloom which caused the main component to fade 
and leave only the baseline noise as desired to enable seeing any small 
signals. The approximate frequency shown in the photograoh was from 
1130 to 1330 MHz. Scan time on Photograph A-2 was 2 seconds/division 
and 20 milliseconds/division for Photograph A-3. Amplitude was 
10 dB/division. 

Photograph A-4 is the spectrum of the transmit pulse with the peak 
approximately 60 dB above the baseline. Some asymmetry is noted but no 
explanation was found to cause this. The spectrum was expected to be 
symmetrical. Photograph A-5 is the same as Photograph A-4 except for 
removal of baseline clipping. Scan time was again 2 seconds and 
20 milliseconds/division. Amplitude was 10 dB/division. 
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Figure A-2. Test Configuration for Obtaining Spectrum Photographs 
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PMOTOGRAPFI A-1. DliTliCTIil) TWIN TRANSMITTIiR 

PULSUS 

SWUF.r 2 MICROSHCONDS/PIVISION 
AMPLITUDE 100 M I M, 1 VOI.TS/Dl V I S I ON 



PHOTOGRAPH A-2. TRANSMIT SPI-XTRUM 

SCAN 1130 MHz to 1330 MHz (20 MHz/DIV.) 
AMPLITUDE 10 dB/DIV. 

SCAN TIME 2 SEC/DI V. 


I 
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PHOTOGRAPH A-3. TRANSMIT SPECTRUM 

BASELINE CLIPPER MINIMUM 

1130 MHz to 1330 MHz (20 MHz/DIV.) 

AMPLITUDE 10 dfl/DlV. 

SCAN TIME 20 mSEC/DlV. 
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PHOTOGRAPH A-4. TRANSMIT SPECTRUM 
BANDWIDTH lOO KHz 
SCAN 0.5 MHz/DTV. 

SCAN TIME 2 SEC/DIV. i 

AMPLITUDE 10 dB/DIV. j 

\ 

i 

1 





PHOTOGRAPH A-5. TRANSMIT SPECTRUM 

BASELINE CL IF PER MINIMUM 
BANDWIDTH 100 KHz 
SCAN 0.5 MHz/DIV. 

SCAN ttmi: 20 mSEC/DIV. 
AMPLI UDE 10 dB/DIV. 
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